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Abstract 

We pursue our study of the antiperiodic dynamical 6-vertex model using Sklyanin’s separation of variables 
approach, allowing in the model new possible global shifts of the dynamical parameter. We show in particular 
that the spectrum and eigenstates of the antiperiodic transfer matrix are completely characterized by a system 
of discrete equations. We prove the existence of different reformulations of this characterization in terms of 
functional equations of Baxter’s type. We notably consider the homogeneous functional T-Q equation which 
is the continuous analog of the aforementioned discrete system and show, in the case of a model with an even 
number of sites, that the complete spectrum and eigenstates of the antiperiodic transfer matrix can equivalently 
be described in terms of a particular class of its (4-solutions, hence leading to a complete system of Bethe equa¬ 
tions. Finally, we compute the form factors of local operators for which we obtain determinant representations 
in finite volume. 


'Univ. de Bourgogne; IMB, UMR 5584 du CNRS, Dijon, France; damien.levybencheton@ens-lyon.fr 
2 ENS Lyon; CNRS; Laboratoire de Physique, UMR 5672, Lyon, France; giuliano.niccoli@ens-lyon.fr 
3 Univ. Paris Sud; CNRS; LPTMS, UMR 8626, Orsay 91405, France; veronique.terras@lptms.u-psud.fr 


1 



Contents 


1 Introduction 

2 The dynamical 6-vertex model 

2.1 The dynamical Yang-Baxter algebra. 

2.2 Left and right representation spaces . 

3 SOV basis in left and right representation spaces 


H 

a 

a 

a 

y 


4 


Diagonalization of commuting antiperiodic transfer matrices 

4.1 The SOV discrete characterization of the spectrum and eigenstates. 

4.2 On the reformulation of the SOV characterization of the spectrum in terms of solutions of func¬ 
tional T-Q equations: homogeneous versus inhomogeneous equations. 


15 

15 


18 


4.3 

4.4 

4.5 


On the reformulation of the SOV characterization of the spectrum in terms of solutions of an 
homogeneous functional T-Q equation: a few preliminary considerations . 

Study of the homogeneous T-Q functional equation: an Ansatz for the Q-solutions. 

On the reformulation of the SOV characterization of the spectrum in terms of solutions of an 
homogeneous functional T-Q equation: proof of the completeness for N even. 


20 

23 


30 


5 Local operators and dynamical Yang-Baxter algebra 

6 Form factors of local operators 

7 Conclusion 

A Theta functions, elliptic polynomials and useful identities 
B Inhomogeneous Baxter equation as reformulation of SOV spectrum 
References 


35 


H 


43 


43 


H 

0 


2 




















1 Introduction 


This paper is a continuation of a previous work of one of the authors |[38l . in which the study of the antiperiodic 
dynamical 6-vertex model was initiated by means of the Separation of Variable (SOV) method introduced by 

Sklyanin i! [501 ED ISO- 

The dynamical 6-vertex model is constructed from an /(-matrix satisfying the dynamical (or modified) Yang- 
Baxter equation ||T9llT5ll . In the context of exactly solvable models of statistical mechanics, it is the archetype 
of interaction-round-faces (IRF) models |5), which describe the interactions of a local variable around faces of 
a two-dimensional square lattice. This is notably the case of the exactly solvable solid-on-solid (SOS) model 
mmmmm, which models the growth of a surface (for instance in the context of a crystal-vapor interface) 
with respect to a flat reference surface. In this context, a height variable is attached to each site of the lattice, 
and the local Boltzmann weights describing the probability of a height configuration around each face of the 
lattice correspond to the non-zero entries of the dynamical /(-matrix. This model plays a crucial role in Baxter’s 
solution Q of the famous eight-vertex model. 

The fact that the /(-matrix of the model depends on an extra parameter, the so-called dynamical parameter 
(related to the height variable in the language of SOS model), results into a modification of the algebraic structure 
underlying integrability compared to what happens in usual vertex models such as the six-vertex model lfl5l 
mm. In particular, as mentioned above, the /(-matrix of the model no longer satisfies the usual quantum 
Yang-Baxter equation, but instead a modified version of this equation, in which the dynamical parameter gets 
shifted by an element of the Cartan. As a consequence, the corresponding Yang-Baxter algebra incorporates 
some additional operator structure acting on the dynamical variable flT8lfT7 1. In practice, the appearance of the 
dynamical shifts may be a problem for actual computations of physical quantities of the model. For instance, 
the partition function of the model with domain wall boundary conditions does not seem to be expressible in 
the form of a single determinant I48ll44l as in the six-vertex case lf22l . As a consequence, the study of the form 
factors and correlation functions of the periodic model in the algebraic Bethe Ansatz (ABA) framework happens 
to be slightly more complicated than in the six-vertex case lf30U3lll32ll . In fact, the latter relies on the use of 
a compact formula, preferably in the form of a single determinant, for the scalar products of Bethe states. The 
problem is that there does not exist, in the ABA framework, a model-independent clear procedure to construct 
(or even guess) such a representation, either for the Bethe states scalar products or for the aforementioned 
partition function, two quantities which are intimately related. 

In the antiperiodic SOS model, which can be solved by SOV llT6l |38l , the situation is somehow simpler. 
First, the space of states of the (finite size) model is finite dimensional for generic crossing parameter, contrary 
to what happens in the periodic case: this means in particular that, when performing a change of basis between 
the canonical basis of the space of states and the eigenstate basis, we only have to deal with finite sums. Second, 
due to the eigenstate representation in terms of separate variables, determinant formulas for the corresponding 
scalar products appear in a much more natural way, somehow intrinsic to the method If38l . We shall see in the 
present paper that such formulas can quite naturally be extended to the form factors. 

The puipose of this paper is three-fold. First, we revisit the study of If38l so as to slightly generalize it 
to cases that correspond to different versions of the SOS solvable model (see [29l l46ll ). related to different 
possible values of the global shift of the dynamical parameter of the model, and that will be useful for our 
further study of the eight-vertex model J4T|. Second, we pursue the study of the spectrum, and discuss the 
reformulation of the discrete SOV characterization in terms of solutions of some homogeneous (respectively 
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inhomogeneous) functional equation of Baxter’s type (and hence of Bethe-type equations), a question that we 
expect to be of primary importance for the consideration of the homogeneous and thermodynamic limits: we 
notably show that, in the case of a model with an even number of sites, the SOV characterization of the spectrum 
is completely equivalent with a description in terms of particular classes of solutions of the usual (homogeneous) 
T-Q equation. Finally, we explicitly write determinant representations for the form factors of local spin and 
height operators. 

The article is organized as follows. In Section[2l we introduce the model and recall the algebraic framework 
for its resolution. In Section^ we construct a SOV basis of the representation space of the model. In Section[4] 
we diagonalize the commuting antiperiodic transfer matrices on the subspace of the representation space which 
corresponds to the actual space of states of the antiperiodic SOS model. We notably discuss the characterization 
of the spectrum and eigenstates in terms of solutions of functional equations. Then, in Section [5] we solve 
the quantum inverse problem for local operators of the model in terms of the antiperiodic monodromy matrix 
elements, and in particular in terms of the antiperiodic transfer matrix. This enables us, in Section | 6 j to write 
explicit determinant representations for the form factors of local spin and local height in the eigenstate basis of 
the transfer matrix. 


2 The dynamical 6-vertex model 


The dynamical 6 -vertex model is associated with a dynamical /?.-matrix of the form 


R(X\t) 


(a{ A) 0 0 

0 e iy,? b(X\t) e iyX c(A|i) 

0 e“* yA c(A| — t) e~ iyv 6 (A| — t) 

\ 0 0 0 


0 

0 

0 


\ 


€ End(Vi <g) V2), 


a(A)/ 


(2.1) 


with Vi — C 2 . Throughout this paper, y G {0,1} is fixed. The R-matrix (12.11) depends on two parameters, a 
spectral parameter A € C and a dynamical parameter t £ to + r/Z, where 1 0 will be specified later. The functions 
a, b, c are given as 


a(X) = 9(X + 77 ), b(X\t) 


9{\)9{t + r 1 ) 

m 


c(A|f) 


9f)9{t + X) 

9 if) 


( 2 . 2 ) 


where 9{ A) = 0\ (A|tu) denotes the usual theta-function (see Appendix lAl) with quasi-periods ir and ttoj ( ( 3oj > 
0). 77 £ C is the crossing parameter of the model, which is supposed to be generic throughout this paper. The 
i?-matrix (12. lb with dynamical parameter t is solution of the quantum dynamical Yang-Baxter equation lfT9lFl5ll 

on V\ ® V 2 <8> Vs, 


Rl, 2 (X 12 \t + 7703 ) -Rl, 3(^13 \t) R 2 , 3 (^ 23 ^ + r t cr 1 ) — R 2 , 3 (^ 23 ^) (A 13 + 770 -I) i?l, 2 (Al 2 |*). (2.3) 

Here and in the following, the indices indicate as usual the space of the tensor product on which the correspond¬ 
ing operator acts non-trivially, and of (o = x, y. z) denotes the usual Pauli matrix acting on Vi ~ C 2 . Also we 
have used the shorthand notation A ij = A, — A j. 

As mentioned in Introduction, the matrix elements of the /i-matrix (12.11) can be understood as the local 
Boltzmann weights of an exactly solvable solid-on-solid (SOS) model (see Figure |T|). In this framework, the 
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dynamical parameter t corresponds to the height variable, and the dynamical Yang-Baxter equation (12.31) is 
simply the star-triangle relation of IRF type 0. The shift (G {±?/}) of the height between two neighboring sites 
can be understood in terms of a spin variable (G {±1}) on the corresponding link. 
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Figure 1: The 6 different local configurations around a face and their associated local statistical weights. 

Remark 2.1. The case y = 1, which was not considered in lf38l . corresponds to a diagonal dynamical gauge 
transformation of the /(-matrix of lf38l of the form: 

R 12 (\ l2 \t) = G 2 {\ 2 \t)G l (\ l \t + i 1 a z 2 ) [Ri 2 (\i 2 \t)] y=0 G 2 (X 2 \t + r 1 a z l )- 1 G 1 (X 1 \t)- 1 , (2.4) 


with 


G(X\t) = e 


( e iy % 

V 0 



(2.5) 


This case is interesting since it enables us to consider a model in which the dynamical parameter is shifted 
by half of the imaginary quasi-period (i.e. to G M + yfw, see Section 12.21 ) so as to recover for instance the 
Boltzmann weights considered in l29l l47ll . It is also useful for the consideration of particular quasi-periodic 
boundary conditions for the 8-vertex model obtained from this SOS model by vertex-IRF transformations Pffll . 

Remark 2.2. It may be interesting to consider the trigonometric limit of (12.1I) - (12.2I) . which corresponds to the 
limit u —>• +?'oc. If to is of the form to = to + yf w, with to independent of u (as considered in this paper, 
see (12.291) 1. then one obtains different limits according to whether y = 0 or y = 1. Up to normalization, the 
trigonometric limit of the case y = 0 corresponds to the trigonometric dynamical 6-vertex model, with 


/ \ \ . s ,/M,x sin A sin(f + rj) sin r] sin(f + A) 

a(A) = sm(A + rj), b(X\t ) =-c(A|i) =-, 

sin t sin L 

whereas, in the case y = 1, one simply recovers the ii-matrix of the usual 6-vertex (or XXZ) modejj). 


( 2 . 6 ) 


4 We use the fact that ie ^ 4 6h(u|cv) — > sin u, and that —>■ e^ v for u, v £ C and 0 < e < 1. 

2 I ' u _> +ioo B 1 (v±enu,\u,) .. ...... ’ 
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2.1 The dynamical Yang-Baxter algebra 


The elliptic quantum group associated with the dynamical /i'-matrix (12.11) was introduced in lfl5l . and its rep¬ 
resentations were studied in lfl 8 l . leading in ifTTI to an algebraic Bethe Ansatz for the corresponding periodic 
SOS model. In this framework, the key object is the monodromy matrix, which provides a representation of 
the corresponding Yang-Baxter algebra. In the case of a tensor product of fundamental representations that we 
consider in this paper, it is defined as the following ordered product of /(-matrices, 

N—1 

M 0 (A|f) = /? 0 ,n(A -£n| t + yy^Ug)--- Ro,i(X ~ £i|t) <E End(V 0 ® Vj <g>... <g> Vn) 

a= 1 

= /A(A|t) B(A|t)\ 

\C(A|f) D(A| t)J [Q] - 

Here N is the size of the model, V n ~ C 2 for n 6 {0,1,..., N}, and £ n , n E {1,..., N}, are inhomogeneity 
parameters. In this context, Vo is usually called the auxiliary space and Vn = ( 8 >„ = , Vj, the quantum space. 
Commutation relations for the entries A(A|t), B(A|f), C(A|f), D(A|t) E End(V|\i) of the monodromy matrix 
(12.71) are given by the following quadratic relation, which is a consequence of (12.31) . 

^o,o'(Aoo '\t + r?S) Mo(Ao|f) Mo'(Ao'|t + V a o) = Mo'(Ao'l^) Mo(Ao|i + V a o') Ro,o'{^oo'\t), (2.8) 

where S is the total 2 -component of the spin: 

N 

s = E<- (2 - 9) 

n =1 

To handle the shifts of the dynamical parameter induced by the relation (12.81) . it is convenient to introduce, as in 
Il38l . some dynamical operators r and T f which commute with local spin operators cr“ and such that 

T±r = (r±r/)T± (2.10) 


This enables us to define a new monodromy matrix incorporating these dynamical operators, 


.Mo(A) = M 0 (A|r) Tt° = 


'A A) B(X)\ 
,C(A) V(X)) m 


( 2 . 11 ) 


where Mo(A|r) coiTesponds to the monodromy matrix (12.71 ) in which we have substituted the dynamical param¬ 
eter t by the operator r, and where 


t^ ct o — 

I T - 


T± 0 
0 Jf 


( 2 . 12 ) 


The operator entries A, B, C, V of the monodromy matrix (12.1 lb act on the space 


®( 6 VD).n = Vn (2-13) 

where D is a representation space of the dynamical operators. Their commutation relations arc now given by 
the quadratic relation 

-Ro,o'(Aoo'l r + ^S) Ado(Ao) Ado'(Ao') = Ado'(Ao') Ado(Ao) -Ro,o , (Aoo , l r )> (2.14) 
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(2.15) 


which follows from (12.81) and from the zero-weight property of the /(-matrix (12.11) : 

Tr aZ ° /?o,o'(A 0 o'|r) T°° T? = R 0 ,o'(X 00 ,\t). 

We also recall the inversion formula for the monodromy matrix: 

Theorem 2.1. The inverse of the monodromy matrix (12.111) is given by the relation 
tV \a~(\ _ ^do „y — r -iynS 


Mo(A) • erg Mq(X - rj ) 0 erg = e" 


■ detqAf(A), 


(2.16) 


9(t + 7/S) q 

in terms of the quantum determinant det q M (A). The latter is a central element of the dynamical Yang-Baxter 
algebra defined as 

det 9 M(A) = a(A) d(A - r/) = e^ s (-4(A) V(X - t?) - 13(A) C (A - t/)) 


= e 


: ” s w - -n) - C(A) B(A - .,)), 


with 


A(A) = lia(A-( n ), D (A) = a(A - 7/). 


(2.17) 


(2.18) 


n=l 


The algebraic Bethe Ansatz ifTTI consists in diagonalizing the transfer matrices of the periodic model, i.e. 
the traces of the monodromy matrix (12.111) . 

T(A) = .A(A) + P(A), (2.19) 

on the subspace of the representation space (12.13b associated with the zero eigenvalue of the spin operator S 
(12.91 ). In this paper, we apply instead the SOV approach to the antiperiodic model (or more generally to a 
K-twistcd antiperiodic model for some parameter n <G C \ {0}), and therefore define the following k - twisted 
antiperiodic monodromy matrices: 

mJ°(A|t) = X^oZ Mo(A|r), A) = M^(A|r)T^, (2.20) 

with X^ = diag(K, kA 1 ). It is easy to see that, thanks to the following properties of the R-matrix (12.11) . 

erf of i?i, 2 (A|f) erf of = /? 1 , 2 (A| —t + ynu>), (2.21) 

x[ K) X^ ] Ri, 2 (X\t) x{ K_1) x[ K ~ 1] = Ri, 2 (X\t), (2.22) 

the monodromy matrices (12.201) satisfy the quadratic relations 


i*o,o'(Aoo'l - r - 7/S + yvru;) M^(A 0 |r) T? l< j (A 0 ,\t) T 7 


o( K ) 


= Mg7(Aor|r)T r °' M^(A 0 |t)T t R 0 , 0 >(X 00 '\t), (2.23) 


-Ro,o'(Aoo' | — r — 7/S + y7roj) Mq K \Xo) Mq^(X 0 ') = My\x 0 ') Mq K \Xo) /?o,o'(Aoo'l r )- (2-24) 

In the following, we will show how to diagonalize, in the SOV framework, the «:-twisted antiperiodic transfer 
matrices, 

r W ( A) = k -1 £>(A) + kC( A), (2.25) 

on some subspace of the representation space (12.131) associated with a particular eigenvalue of the operator 
S r = 7/S + 2 r. (2.26) 
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2.2 Left and right representation spaces 


In the following, we may particularize by a subscript L (respectively 1Z) the left (respectively the right) represen¬ 
tation spaces for the spin and dynamical operators. Hence, for instance, the operator entries of the monodromy 
matrix (12.111) acj^to the right on the “ket” space 1 C/^ VD) n = B( 6V d),n an d to the left on its restricted dual space, 
the “bra” space ID^ VD) N - 

j s an i n (j n it c dimensional representation space of the dynamical operators algebra (12.101) with r- 
eigenbasis left (covectors) and right (vectors) respectively defined as 

(t(a)\ = <t(0)|T~“, \t(a)) = T“|f(0)), Va G Z, (2.27) 

such that 

(t(a)\r = t(a)(t(a)\, r|f(a)) = f(a)|f(a)), t(a) =—rja + to, Va £ Z, (2.28) 

with the normalization (t(a)\t (&)) = 6 a ,b, Va, 6 G Z. In this paper, we fix the value of to to be 

*0 = ~2 N +X f +y ? W ’ With x € {0,1}, 

such that (x, y) / (0, 0) when N is even. We denote the corresponding representation space as 

cnz 

We respectively denote the left and right spin basis in V n as 

(n, h n \a z n = (1 - 2h n )(n,h n \, cr z n \n, h n ) = (1 - 2h n )\n, h n ), h n € { 0,1}, (2.30) 

with (n, h n \n, h! n ) = Sh ri .,h'„ f° r an y n € {1,..., N}. Hence the states 

(®n=i( n i h n \) <S> (t(a)|, (®n = i\n, h n )) ® |t(a)), (2.31) 

obtained by tensoring common eigenstates of the commuting operators r and a*, 1 < n < N, provide left 
and right dynamical-spin basis in D^ VD i n and D^ vd . N , respectively. The following scalar product is thereby 
naturally induced in the linear space D’^ VD j n : 

N 

( \n,K) ® \t{ a )),®n=l\ n , h n) ® l^a'))) = 8a t a> JJ 5 hn ,h' n - (2-32) 

72—1 

Remark 2.3. In the case (x, y) = (0, 0) in (12.291) . the corresponding representation of the dynamical Yang-Baxter 
algebra is well defined only when N is odd (to avoid singularities in (12.21) ): this case was analyzed by SOV in the 
paper 11381 (see also lf!6l ). The analysis that we present here will address both the case N odd for (x, y) = (0,0) 
and the case N even and odd for (x, y) ^ (0,0). This enables us to consider more general SOS models (such 
as for instance as in Il29ll46l ). and will be useful for our study of the eight-vertex model with different types of 
quasi-periodic boundary conditions lf4TI . 

5 These notations may be confusing since ®(6VD), N corresponds in fact to a left-module (and D^ V d) n to a right-module) of the elliptic 
quantum group, but they agree with some more physical convention also used in 1381 . 


(2.29) 


(x,y),N 









(2.33) 


The operator S r (12.261) defines a natural grading on Tn ’ £ ^" R 


(6VD),N' 


^(6VD),N ttJ r=-oo- lu; (6VD),N’ 


where ^(gvD^N LS the 2 N -dimensional linear eigenspace coiTesponding to the eigenvalue 2rrj + xtt + ynu of S r . 

In terms of the dynamical-spin basis the linear (covector) space D| 6VD ; n and the linear (vector) space ©| 6VD ) N 
are respectively generated by the elements 

(<8>^ =1 (n, h n |) (8) (f r , h | with (t rjh | = (*o,h|T^ , (2.34) 

and 

((g>^ =1 | n, h n )) (8) |f r , h ) with |f r ,h) = Tf r \t 0 ,h)- (2.35) 

Here we have set 


U,h = -|s h + x^ + y^w + rr), with s h = ^(1 - 2 h k ) and h = (hi,... , h N ). 


7T 7T 

;-h y— 

2 2 


(2.36) 


fc=i 


• • r 77 — /'Rs) 

Proposition 2.1. For each r G Z, the finite-dimensional vector spaces /, are invariant under the action 


of the operators 

A(A|t), D(A|r), B(A), C(A). 


(6VD),N 


(2.37) 


Proof. The commutation relations 

[A(A|t), S] = [A(A|r),r] = [D(A|r), S] = [D(A|t),t] = 0, (2.38) 

[B(A), S] = 225(A), [B{ A), r] = - V B( A), (2.39) 

[C(A), S] = — 2C(A), [C(A), r] = V C( A), (2.40) 

imply that 

[A(A|r), S T ] = [D(A|t ), S T ] = [B(A), S T ] = [C(A), S T ] = 0, (2.41) 

which means that ^ arc invariant under the action of these operators. □ 


For a SOS model with free boundary conditions, the physical space of states corresponds to the whole rep¬ 
resentation space B(6 vd).n- ^ we impose different types of boundary conditions, the physical space of states 
will correspond only to a subspace of B( 6 V d),n- F° r instance, in the case of periodic boundary conditions, it 
corresponds to the subspace of (12.131) associated with the zero eigenvalue of S, whereas in the case of antiperi- 
odic boundary conditions that we consider in this paper, it corresponds to the subspace of (12.131) associated with 
the eigenvalue X7r + yTrcc of S r , i.e. to D^ d n . In fact, we shall see in Section [4] that the commutation of the 
antiperiodic transfer matrices is ensured on this subspace only. Nevertheless, the construction of the SOV basis 
that we present in the next section holds in the whole representation space B( 6V D),i\i- 

Remark 2.4. In the case of periodic boundary conditions, the space of states of the (finite-size) model is usually 
infinite dimensional, which may be a technical inconvenience for the study of the model. In fact, due to this 
reason, one usually deals with restricted models for which the crossing parameter is rational and the space of 
states is finite dimensional (such as for instance in the original paper [0). Let us stress that, for antiperiodic 
boundary conditions, we do not have this problem since the space of states D^ d . N has dimension 2 N . 
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3 SOY basis in left and right representation spaces 


For usual vertex models such as the 6-vertex model Il50ll5ni43ll2lll39ll40ll42ll . the SOV approach to diagonalize 
the transfer matrix is based on the construction of a basis of the space of states of the model which explicitly 
diagonalizes the action of one particular generator of the Yang-Baxter algebra. In the dynamical case, however, 
the corresponding SOV basis for the antiperiodic transfer matrix fl6ll38l only partially diagonalizes the operator 
T>(\) (or A(X)). In other words, the latter still acts as a shift operator on some “dynamical” part of the corre¬ 
sponding basis, whereas the operators B( A) and C( A) act as a sum of shift operators on the “spin” part of the 
basis. In this section, we present the explicit construction of this SOV basis in the left and right representation 
spaces of the model. 

In each of the subspaces B^ 6VDj N , we define the following reference states : 



(3.1) 


where we have used the notations 0 = (hi = 0,... , =0) and 1 = (hi = 1,..., = 1), and where N is a 

fixed normalization constant that will be specified latter. It is easy to see that these references states are D(A|r) 
and A(A|r)-eigenstates with eigenvalues given in terms of a(A) and d(A) (see (12.181) ). 



(r, 0| A(A|t) = a(A) (r, 0|, (r, 0| D(A|r) = e 


(3.2) 


D(A|t) |l,r) = a(A) |1 ,r), A(A|r) |l,r) = e 


(3.3) 


and are annihilated by the action of the operators £>(A). Then, for each N-tuple h = (hi ,..., /in) €= {0,1} N , 


we construct a state (r, h| € ID>^ VD) n an ^ a state |h> r) € 0^ VD j n as 



(3.4) 



(3.5) 


Remark 3.1. It is easy to see that, Vh 6 {0,1} N , 


(r, h| r = t r>h (r, h|, r |h, r) = t r>h |h, r), 

(r, h| S = s h (r, h|, S |h, r) = s h |h, r), 


(3.6) 

(3.7) 


with f r h and Sh given by (12.361) . which also implies that (r, h| € ^(6 vd) n anc ^ l^ 1 - r ) e ^fevo') n 


We have the following result: 

Theorem 3.1. Let I = 7rZ + 7rwZ and let the inhomogeneity parameters £i,..., £n £ C be such that 


VeG {-1,0,1}, ^a-^b + ep^T if a f b. 


(3.8) 
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Then the set of states {(r, h|, h £ {0,1} N } defined by (13.41) form a basis o/B^vd) i\t Similarly, the set of states 
{|h, r),h £ {0,1} N } defined by (13.51 ) form a basis of N . Moreover, the action ofV(X) on the states 

( r , h| in B^ VD) N (respectively on the states |h, r) in B^ VD ^ n ) is given by 


(r, h| V{ A) = d r _i ih (A) (r - 1, h|, 
T>(\) |h, r) = d r+ i ih (A) |h, r + 1), 


where 


dj-,h(A) — e 


6{tr, h) 


Q{t r ,i 


with ^) = ^ a - r] h a 


n =1 


(3.9) 

(3.10) 


(3.11) 


Proof. The action of the operator 'D(X) on the states (r, h| and |h, r) can respectively be computed by means 
of the following dynamical 6-vertex commutation relations issued from (12.141 ): 


C(p)V( A) 


V(\)C(n)0(\-ti + ri)0(T) 


- V(n) C{ A) e iy(A -' i) @(r ? )0(r + A - /r) 


e iyv 

9( X — h)6(t 


vV 


(3.12) 


V(X)C^) 


e -iyv 

6(n - X)6 (t + r]) 


0{n - A + r})0{T)C{n) V(X) 


- e -*y(^-A) e ^)0 ( T _ M + A) C(A) V(ff) 


(3.13) 


Hence, using the fact that d(£„) = A (f n — rj) = 0, n = 1,..., N, so that only the direct terms contribute in the 
above relations, we obtain (13.9b and (13.101) . 

It remains to prove that the states (r, h|, h £ {0,1} N , form a set of 2 N independent states, i.e. a basis of 
fi(6VD) N (the fact that the states |h, r), h £ {0,1} N , form a basis of N can be proven similarly). Let us 

suppose that we have a relation of the form 

^2 c h (r,h|=0, (3.14) 

he{o,t} N 

for some set of complex coefficients Ch £ C, h € {0,1} N . Let k = (k\,..., k n) be any given N-tuple in 
{0,1} N . Then, by applying the product 

N 

with k n £ {0,1}, k n = k n +1 mod 2, (3.15) 

72=1 

to (I3.14h . we get the identity 

N 

C k J] d r -nM{f ( t y ) {r - N,k| = 0. (3.16) 

72=1 

Since d r _ n k(^n' r ^) 7 ^ 0, 1 < n < N, and (r — N, k| f 0, this implies that c k = 0, which ends the proof of 
Theorem l3.ll □ 
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We want to determine the action of the remaining generators of the Yang-Baxter algebra on the basis el¬ 
ements (13.41) and (13.51) . and in particular of the operators -8(A) and 8(A) constituting the antiperiodic transfer 
matrix. To do this, we shall rely on the following lemma: 

Lemma 3.1. V( A), -8( A) = e~ iyX B( A) and C( A) = e* yA C(A) are entire functions of X which satisfy the follow¬ 
ing quasi-periodicity properties: 


V(X + tt) 

= (-1) n P(A), 

V(X + 7T w) = 

e 2 *E^=i(?n-J?/2)+ir;S J)(X), 

(3.17) 

B(X + n) 

= (—1) N 8(A), 

B( X + TTW) = 

J g 2 ’En=iK» _ W 2 ) _ *s T -8(A), 

(3.18) 

8(A + 7r) 

= (-1) n C(A), 

8( X + nw) = (~e~ 

| N e 2i T,n=l(^-vfy)+iS-r C(X). 

(3.19) 


Proof The fact that T>{ A), 8(A) and C(X) are entire functions of A is a simple consequence of the definition 
of the dynamical 6-vertex monodromy matrix in terms of the dynamical 6-vertex ii-matrix (12.11) . The quasi¬ 
periodicity properties (13.171) of T>( A) follows directly from the explicit form (I3.9I ). (13.101) . (13.1 II) of its action 
on the basis (13.41) or (13.51) and from the quasi-periodicity properties of the usual theta function 9( A) = 0\ (A|a>). 
The quasi-periodicity properties (13.191 ) of 8(A) can be deduced from the quasi-periodicity properties (13.171) of 
V(X) by comparing the quasi-periodicity properties in A of the two members of the commutation relation (13.121 ) 
rewritten in the form 


8( A) V{p) = e -iy(A-A*-u) 


<9(A-/x)<9(r + r?) 
9(X- p + r)9(q) 


V(X)C(p) 


9(X - p - ij) 9 (t) 
9{X- p + r)9{q) 


C(p)V(X). 


Similarly, the quasi-periodicity properties of -8(A) can be obtained from (13.171) and from the commutation 
relation 


-8(A) V{p) = 


9(X - p + q)9{r + q S) 
9(t + qS - X + p) 9(q) 


B{p)V{ A)-^-^ 


9(X - p) 9 {t + qS - ij) 
9 (t + qS - X + p) 9(q) 


V(X)B(p). 


□ 


Theorem 3.2. Under the hypothesis (13.81) of Theorem 13.71 the action of the operators A(A|r), D(A|r), -8(A) 
and C(X) on the basis elements (r, h| g/B^^) N (respectively on the basis elements |h, r) o/B|g(^ N ) take 
the following form: 

I) Left representations: 


(r,h| D(A|r) = d,._i ih (A) (r,0| J][ 


N / T-C(£ n )T+^" 


n= 1 


D(Cn - V) 


(3.20) 


(, h| C(A) = ± IJ T >l. (3.2D 

(r . h| B (A) = £ n g g;_ d g, ) 4U(d‘-'- ) ) T„-h|, <3.22, 


a= 1 


and the action o/A(A|r) is completely determined by the quantum determinant relations (12.171) . 
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II) Right representations: 


D(A|t + rj) |h, r) = d r+1)h (A) 


’it /^ T+C^n-r?) 

,Mv °&~v) ) 


r ) 


C (A) ih,r> - £n '<^2, “(e-») r.-h, r > 


a=l 

N 


0(Zr,h) 

e( h a)_^\ 9(t _ r _h — A + £a °^) 


b+a J 


m IM - £ ' n 2(221 VA.hKM |T>, r> 

a= 1 


#(/-r,h) 


Q(td l a) c(h-b)\ 

b+a ~ J 


(3.23) 

(3.24) 

(3.25) 


and the action of A(A|r) /.v uniquely determined by the quantum determinant relation (12.171) . 

In the above expressions, we have used the definition (13.111) of d r h(A) and ^i a \ and we have set 
Tq (h \,..., /in) = (hi, • • •, h a ± 1 ,..., /in), 


-£!,,hW = (-1)*+^ e 2iyr i fffff a(A). 


6»(/_ rjh ± r/) 


(3.26) 

(3.27) 


Proof. The action 17.201) . (13.231) of D(A|r) = D(A)T+ and D(A|r + rj) = JfV(X) on (r, h| and |h, r) 
respectively follow directly from the formulae (13.91) . (13.101 ) and from the commutation relations (12.101 ). 

To determine the action of £>(A) and C(X) on a left state (r, h| (respectively on a right state |h, r)) associated 
with a given N-tuple h = (hi, ...,/in), we first compute these actions at the N special points A = HI 
n = 1,..., N. From the definitions (13.41 ). (13.51) . and the fact that the product C(f n ) C(f n — rj) vanishes when 
evaluated at one of the inhomogeneity parameters f n (see (15.91 )1. we get 

(r,h\C(f^) = d^-M) <r,T+h|, C(H" } ) |h ,r) = d(H°) |T“h,r). (3.28) 

The action of on the left and right states (13.41) and (13.51) can be computed by means of the dynamical 

Yang-Baxter commutation relations (12.141) and the quantum determinant relations (12.171) . We obtain 

<r, h| = -Affi 1 -' 1 "!) {r,T;h| 


and 


= (-l)* +y+ * y <»”> ffx f (r,T-h| 

V(t-r ,h - 1)) 

B((F) |h,r> = -A«<‘»>) | T +h, r) , 

= (-1)«»+^ A(d k ">) e 2 ‘y'” / (t '' h |T+h, r), 

u(Z-r,h + V) 


(3.29) 


(3.30) 


in which we have used the quasi-periodicity properties of the theta function. 

Finally, the action of B( A) and C(A) to the left 13.2 1 1 - 13.221 ) and to the right (I3.24I) - (I3.25I) are obtained from 
(13.291 ) and (13.301) and from Lemma [Til by means of interpolation formulae (see (IA.6I0 . □ 

Remark 3.2. It is easy to check that the representations of the Yang-Baxter generators in the SOV basis coincide 
with their representations in the so-called F-basis If33l as obtained in |[Q. Such a relation between the SOV 
construction and the F -basis was first noticed in 
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The basis elements (13.41) and (13.51) have been constructed independently in the left and right representation 
spaces respectively. It is nevertheless not difficult to compute their scalar product f38ll : 

Proposition 3.1. Under the hypothesis (13.81) ofTheorem \3.1\ let us define, for each h E {0,1} N and r E 7L, the 
N x N matrices 0' r,h ' 1 of elements 


[ 0(r,h) ] ij = 0J-l(£i Ai) _ £r), with £ r = ^ ^ & + tr, 


where the functions t), are defined as 


•dj(X) = ^2 e i7rNw(n+ 2“N) + 2iN ( n +2-M)( A -D ) 0 < j < N - 1. 

7tEZ 

Then, one can fix the normalization constant in (13.11) such that, Vh, h E {0,1} N , Vr, r E Z, 

0(t r ,i) e - iy ^=i hfe 


(3.31) 


(3.32) 


(r, h|h, r) = S r ^ <5] 


h.h 


det|\| [0( r . h )] 


(3.33) 


Proof The fact that (r, h| and |h, f) are both eigenstates of S r with respective eigenvalues 2 rrj + X7r + yiru; 
and 2 frj + x-zr + yiru implies that (r, h|h, f) vanishes if r f f. Similarly, from the consideration of the matrix 
element 

(r + 1, h| T>(\) |h, r) (3.34) 

and the explicit expression (13.91) . (13.101) and (13.1 II) of the action of 'D(X) to the left and to the right, one gets that 
(r, h|h, f) vanishes if h h. One also gets that 

(r + 1, h|h, r + 1) _ 0(i r , h ) 9(t r + 1 , 1 ) 


(r, h|h, r) 9{t r ,i) 0(t r +i,h) ’ 

Also, from the consideration of the matrix element 

<r,h|C(£ 0 )|T+h,r) 


(3.35) 


(3.36) 


for any N-tuple h such that h a = 0, and the explicit expressions (13.281) . (13.241) of the action of C(fi a ) to the left 
and to the right, one gets that 


>T+h|T+h,r) 
(r, h|h, r) 


0(t r , h ) jj0(d O) -^ 6) ) 




(3.37) 


Since {??j_i}i< 3 <N is a basis of the space of theta functions of order N and of norm 0 (see Appendix lAl). we 
can use (1A.7I) to express the determinant of the matrix 0 !r ' h l as 


det [0( f ’ h )] = c N 0(fr,h) Ue^ -$ hi) ), 


(3.38) 


i<j 


where cm is a constant (i.e. it does not depend on the variables on r nor on h). Hence, it follows from (13.351) . 
(13.371) that 


(r,h|h,r) = c N 




det N [0( r - h )] 


where cn is a constant. 


(3.39) 

□ 
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It follows from Proposition 13. 1 I that 


_ ^ |h, r)(r, h| 

(r, hlh, r) ’ 

he{o,i} N 


(3.40) 


with normalization (r, h|h, r) given by (13.331) . provides a decomposition of the identity on the whole represen¬ 
tation space Ib^gvDj i\r Similarly, for any given r£ Z, 


I r 


£ 

he{o,i} N 


|la, r)(r, h| 
(r, h|h, r) 


(3.41) 


provides a decomposition of the identity on the subspace B 


(r,C/K) 

(6VD),N 


of B 


C/TZ 

(6VD).I\T 


4 Diagonalization of commuting antiperiodic transfer matrices 


In this section, we diagonalize the K-twistcd antiperiodic transfer matrices (12.251) in the space of states 
of the antiperiodic model. As usual in the SOV framework, the spectrum and eigenstates are completely char¬ 
acterized by a discrete system of equations involving the inhomogeneity parameters of the model. We explain 
how to rewrite this discrete system in terms of functional T-Q equations of Baxter’s type, so as to obtain a new 
characterization of the spectrum and eigenstates in terms of solutions of Bethe-type equations. 


4.1 The SOV discrete characterization of the spectrum and eigenstates 

The diagonalization of the antiperiodic transfer matrices in B^^^ has been performed in Ifl6ll38l in the case 
(x, y) = (0,0) with N odd. The corresponding SOV procedure can easily be extended to the more general cases 
that we consider here, and we have the following result: 


_ 

Theorem 4.1. For any fixed k £ C\ {0}, the n-twisted antiperiodic dynamical 6-vertex transfer matrix T (A) 
(12.251) defines a one-parameter family of commuting operators on All these families are isospectral, 

i.e. the set of the eigenvalues ofT^\ A) is the same for all the values of k G C \ {0} and we can denote it 
with £~=. Moreover, for any fixed H-tuple of inhomogeneities (£i, ■ ■ ■, £n) £ satisfying (13.81) . the spectrum 
ofT (h} (A) in is simple and coincides with the set of functions of the form 


INI 

t(A) = JV y (*“- A > 
< 2=1 


flQo.O ~ A + jq) -r-r 0 (A — &,) - . 

m,o) UaJ ’ 


(t(6 ),-..,t(^)) eC N , 


(4.1) 


which satisfy the discrete system of equations 


t(Ca)t(& -v) = (-l) x+y+xy A(^)D(^ a - qf Va € {1,..., N}. 


(4.2) 


The right T^\x)-eigenstate |\D^) 6 B^^^ N and the left T^(A )-eigenstate (T^l G ^(evD) N associated 
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with the eigenvalue t(A) G Y^cire respectively given by 


i*r’>= E n 

hG{0,l} N o=l . 


K ~ 1 e i y r i 




(h a ) 


<i , t“ ) i= e n[(—)N: 

hG{0,l} N o=l 


Dfe - r)) 

det [0(°' h )] (0,h|, 


det r©(° h )1 |h, 0), 
N L J 


(h a ) 

a 


where the coefficients q[ A,y are (up to an overall normalization) characterized by 


_ D {ia-rf) _ t(£a) 


% 


( 0) A x,y(?a) : 


with A x , y (A) = (-l) x+y+xy A(A). 


(4.3) 

(4.4) 


(4.5) 


Proof. Let us first recall that, from Proposition 12.11 the action of the twisted antiperiodic transfer matrix t' 1 ' 1 (A) 

(6VD),N ^(6VD),N- II follows from (EH that 


preserves the subspace of 


T W ( Ai) r (K) ( A 2 ) = tr 12 \r 1i 2 (X 12 \-T-r, S + yvrtu)” 1 A4^(A 2 ) ^^(AO i?i, 2 (Ai 2 |r) 


= tr i2 

= tr 12 


«i, 2 (Ai 2 | -r-rj S + yttw)- 1 M^A^r + r?af) M^(Ai|r) T^T" 2 i?i, 2 (Ai 2 |r) 
i?t ) 2 (Ai 2 |r) J R 1 , 2 (A 12 | - r - pS + yTTu J )- 1 M^\X 2 )M[ K) (Xi) 


■(«)/ 


(4-6) 


where we have used both the cyclic property of the trace and the zero-weight property (12.151) of the /^-matrix. 
Hence, (14.61) acts as 


tr i2 


M K) (A 2 )mS k) (Ai) 1 = r W (A 2 )T W (A 1 ) 


(4.7) 


on the left subspace B^y^ n of B^ VD ) N associated with the eigenvalue xtt + yiruj of S T = pS + 2r. Com¬ 
mutativity on the right subspace B^^*. N of B^ VD) N follows by inserting the decomposition of the identity Iq 

(13.411) . Hence, for any fixed k € C, the r,'-twisted transfer matrices T ' (A) define a one-parameter family of 
commuting operators on B. 

It follows from the quasi-periodicity properties of the operators B( A) and C( A) (Lcmma l3. 1 1) that the restric¬ 
tion T^ K \X) of e* yA T^\x) on is a theta function of order N and of norm a\ = ^ + fo,o (see 

Appendix lAl). which means that the action of T^\ A) on any vector of ' s completely determined by its 

action at N independent points with respect to at (see ( IA.61) ). It also means that any eigenvalue function r K) (A) 
of T (k) (A) on b|qvd^ takes the form 


N 

t w (A) = ^yy«L-A) 

a— 1 


^(fp,0 - A + £a) TT &(X - ib) \ 

*(*>, o) ll 9(i a - 6) Uaj; 


(4-8) 


in terms of some N-tuple of complex numbers (f( K )(£i),..., t( Ac ^ (Cn)) - The condition for a function t^ K ^(A) of 
the form (14.81 ) to be an eigenvalue of T^\X) on B^^^ N is equivalent to the fact that there exists a non-zero 
vector 



E 

he{o,i} N 


^ K) (h) 


(0,h| 

(0, h|h, 0) 


G B 


(0 ,C) 

(6VD),N ’ 


(4-9) 
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with i/>^(h) = |h, 0), such that 


(4.10) 


Vh € {0,1} N , (^ k) |T (k) (A) |h,0) =t( K) (A)(^ K) |h,0). 

By computing the action of (A) on |h, 0) at the N independent points £n' r * by means of (13.281) . (13.301) . one 
obtains that the condition (14.101) is equivalent to the system of equations 


*(") (h ) = (_i)x+y+xy a(^)) ®W(T+h) + kd(^>) ^ K) (T~h), 


(4.11) 


for any n G {1,... , N} and h G {0,1} N . Taking into account the fact that A (Fn *) = l)(if*, 0 ' 1 ) = 0, one can 
rewrite this system of equations as the following system of homogeneous equations: 


(K) id 0) ) -K - 1 (-i) x + y +* y a( d 0) )) f #(h)\ = /o 


-k D(d^) 




#(T+h )J V°y 


(4.12) 


for any n G {1,... , N } and any h G {0,1} N such that h n = 0. It follows that a function t ( - K> (A) of the form 
(14.81 ) is an eigenvalue of T^\ A) on N if and only if this system admits a non-zero solution, i.e. if and 

only if the determinants of all 2 x 2 matrices in (14.121) vanish: 


(4.13) 
can be written 


t w ({< 0) )tMgM) = (-l) x+y+ ’‘ y A(4°>)D(4 1 )), Vn€ {1.....N}. 

The condition for a function t^(A) of the form J4.KI to be an eigenvalue nl T ' (A) on ILdff 1 N 
similarly and one obtains the same system of equations (14.131) . The announced isospectrality is then a trivial 
consequence of the fact that the conditions (14.81) and (14.13b are the same for all the values of k G C\{0}. 
Hence, we can omit the upper indices (k) and denote the eigenvalue functions of T (A) simply by t(A). 

Finally, it is easy to see that a given function t(A) of the form (14.81) satisfying (14.131) corresponds to a 
unique eigenvector (\k^| of T^\ A), so that the spectrum of T°' * (A) in N is simple. Indeed, since 

A (?n°^); D (£n^) X 0, the solution of (14.121) is uniquely fixed (up to an overall normalization) by the requirement 

#(T+h) 


V>t (K) ( h ) 


= (|_x) x+y+xy 


A(^ 0) ) 


(4.14) 


for any n G {1,... , N} and any h G {0,1} N such that h n = 0. In other words, it means that this eigenvector is 
given by the factorized formula (14.4b - (14.5b . The proof for the right eigenstates is similar. □ 


As usual within the SOV approach, the eigenstates of the transfer matrix are obtained as separate states 
on the SOV basis, so that their scalar product can straightforwardly be expressed as a determinant issued from 
(1X331) : 


= det fop] , 

where J-ip denotes the N x N matrix of elements 




h=0 


A x,y(£q) 

D(£a - 77) 




(4.15) 


(4.16) 


Such a representation can be obtained for any t(A), t'(A) G Note that it does not depend on the value of the 
twist K. 
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4.2 On the reformulation of the SOV characterization of the spectrum in terms of solutions of 
functional T-Q equations: homogeneous versus inhomogeneous equations 

_ f K \ 

It follows from the previous study that the spectrum of T (A) is given as the set of all entire functions t(A) 
satisfying the quasi-periodicity properties 

t(A + vr) = (-l) N+y t(A), (4.17) 

t(A + TTU) = (_ e -2iA-i 7 ro J )N e 2i[ELl^-^+ x f]t(A), (4.18) 


and such that, for each n G {1,..., N}, the matrix 


= 1,-DKi 1 ’) t««) ) 

is of rank one, i.e. that there exists a non-zero vector 


Qt,n 


(4.19) 



(4.20) 


which satisfies Di n ■ n = 0 and which can be used to construct the corresponding eigenstates. In other words, 
the system of quadratic equations (14.21) which completely characterizes, together with the functional form (14.11) 
of the eigenvalues, the transfer matrix spectrum S^=, is equivalent to the following condition: 


Vn G {1,... ,N}, 3 (q?°i>qS) i 1 (0, o) such that 

Vh n G {0,1}, t(^) qgg = A x ,„(&*">) qg n " +1) + d(^) qg^. (4.21) 

Hence the system of equations (14.211) coiTesponds to a discrete version of Baxter’s famous functional T-Q 
equation f6j. However, in its present form, this characterization of the spectrum in terms of a discrete set of 
equations which strongly depends on the inhomogeneities of the system (subject to the condition (13.81) ) does 
not seem very convenient for the study of physical quantities of the model, since it does not allow for an easy 
determination of the homogeneous and thermodynamic limits. On the contrary, Baxter’s T-Q equation in its 
usual (i.e. functional) form is smooth with respect to the homogeneous limit and, thanks to the equivalence with 
a system of Bethe equations, makes it possible to use some standard techniques to study the thermodynamic 
properties of the model under consideration. It is therefore important to be able to pass from the discrete to the 
continuous picture or, in other words, to find an equivalent reformulation of the SOV discrete characterization of 
the transfer matrix spectrum and eigenstates in terms of some particular class of solutions on the whole complex 
plane C of a functional T-Q equation of Baxter’s type. Note that the existence of such a reformulation has been 
already proven for several integrable quantum models solved by SOV 1351136 ,, .211, notably for the antiperiodic 
XXZ spin chain 1421 which constitutes a limiting case of the present model (see Remark [2~2l i. 

Hence, the problem one wants to solve can be formulated as follows: does it exist, for each t(A) G a 
function Q( A) on C, in a class of analytic functions that has to be precisely determined, such that t(A) and Q( A) 
satisfy the continuous (functional) version of (14.211) : 


t(A) Q(A) = a x , v (A) Q (A - ??) + d(A) Q (A + v ). 


(4.22) 
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If moreover, for each n £ { 1 ,... , N}, (Q(£n),Q(£n ~ v)) 7 ^ (0, 0 ), then this solution Q(\) provides, through 
the identification q^ 1 ") = Q l " 1 ), all the vectors (14.201) satisfying the condition (14.211) and enabling ones to 
construct the corresponding eigenstates through (I4.4I )- (I4.3I ). Provided Q(X) can be factorized in a generic form, 
the determination of all eigenvalues t(A) £ s r can therefore be reduced to the determination of the set of roots 
of the corresponding solution Q( A), i.e. of the solution of a system of Bethe-type equations. 

For any given t(A) £ E^=, the equation (14.221 ) is a second-order finite-difference equation which may in 
principle admit two independent solutions Q( A). The whole problem is therefore to determine what could be 
the functional form of these solutions, and whether this form does or not depend on the particular t(A) £ S T 
we consider (problem of the completeness of the associated system of Bethe equations). In general, this may 
not be an easy task, since the functional form of the Q-solutions to this finite-difference equation may be quite 
different from the functional form of its coefficients: for instance, in the present case, it is obvious that (14.221 ) 
cannot admit, for t(A) £ Ef, any solution of the type 

M 

Q(A) =c Q e aX n 0 (A-Ay), cq± 0 , M € N, a, Ai,..., A M G C, (4.23) 

3 = 1 

since the terms t(A) Q(A), a(A) Q(\—rf) and d(A) Q{X+rj) would all have different quasi-periodicity properties 
and would therefore be linearly independent. A possible (and quite usual) way to solve the problem would be to 
explicitly construct the so-called Q-operator E1H1 (see also for instance 1451 17. 3ll8l fT2l ) and to determine the 
functional form of its eigenvalues. This procedure may however be quite involved and, to our knowledge, has 
never been performed in the case of the antiperiodic dynamical 6 -vertex model. 

It was recently suggested in the context of the so-called “off-diagonal Bethe Ansatz” iflOl [9j] that one may 
avoid these difficulties by considering, instead of (14.221) . a generalized functional equation. The idea is to allow 
some freedom in the rewriting of the discrete equations into the continuous one, and in particular the presence 
of an inhomogeneous (“off-diagonal”) term, so as to force the latter to admit solutions of the form (14.231) . 
The equivalence of the discrete SOV characterization of the spectrum and the solutions of such generalized 
functional T-Q equation is then quite simple to prove, and it is probably the easiest way, in the context of SOV, 
to obtain a complete system of Bethe-type equations (see for instance If24ll421 ). Concretely, in the present case, 
such a reformulation would rely on two main ideas. The first one is that one can in fact allow, in the rewriting 
of (14.211 ) into some functional finite-difference equation, some modification of the coefficients a(A) and d(A) 
by a gauge transformation of the form 

A(A) = /(A) a(A), d(A) = (4-24) 

which leaves the product a(A) d(A — 77 ) (and therefore the set of equations (14.21) 1 unchanged. One can therefore 
choose the function /(A) such that, for functions Q( A) of the form (14.231) . the three terms t(A) Q( A), A(A) Q(X— 
77 ) and 13(A) Q (A + 77 ) obey the same quasi-periodicity properties with respect to the periods 7 r and nco. The 
second idea is that, in the rewriting of the discrete finite-difference equations into a continuous (functional) 
one, one can in fact allow the presence of an extra term as long as the later vanishes at each of the points ^n \ 
n £ {1,..., N}, h £ {0,1}, so that the discretized version of this equation effectively coincides with a gauge 
transformed valiant of (14.211) . The resulting equation will then be of the form 

t(A) Q( A) = /(A) A Xjy (A) Q (A - 77 ) + Q (A + 77 ) - a(A)d(A)F(A). (4.25) 
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In fact, there exist many possible rewriting of this form, depending on the function /(A) and the degree M of 
the class of solutions Q(X). A consistent choice for /(A) and Q(X) is for example: 


/(A) = /M(a)= r v** 

M 0(X- n + tofl) 


M 


Q(A) = n^(A-A i ), 

3= 1 


(4.26) 


in terms of some arbitrary complex parameters ft and // and of some roots Ai,..., Am- With such a choice with 
fixed parameters ft and //, it is natural to expect that the reformulation of the transfer matrix spectrum in terms 
of the solutions Q( A) of degree M of (14.251 ) is complete as soon as M > N 0 . In the appendix iBl we explicitly 
detail the case M = N, whose result can be summarized as follows: 


if M = N, for anjO fixed // <G C and ft <G C\ {0}, then t(A) € if and only if there exists a function Q(X) 
of the form (14.261) . solution with t(A) of the inhomogeneous functional equation (14.25b for /(A) = (A) 

given by (14.26b , such that (Q(ft n )-, Q{£,n — v)) 7^ (0> 0) f°r eac h ^ € {1, • ■ • , N} (see Theorem IB. lb . 

The function F( A) = F^q(X) in (14.25b is completely determined in terms of /(A) = fjf\ A) and Q( A) by 
imposing that the r.h.s. of (14.25b is an elliptic polynomial. 

Such a functional equation reformulation of the transfer matrix spectrum has the clear advantage to be 
completely smooth with respect to the homogeneous limit in which all parameters ft n tend to the same value 
(which was obviously not the case for the initial formulation). It moreover enables us to rewrite the transfer 
matrix eigenstates obtained by SOV in a form very similar to a Bethe vector, i.e. as some multiple action of 
the operators V, evaluated at the roots of the Q-function (the “Bethe” roots), on some specified pseudo-vacuum 
state (see Appendix |Bj Corollary IB. lb . Nevertheless, the possibility to use it as an efficient tool to analyze also 
the thermodynamic limit is still an open question. Indeed, there is an important price to pay for the apparent 
simplicity in the deriving of the functional equation: the presence of the inhomogeneous term F( A) makes the 
resulting Bethe-type equations a priori much more difficult to solve than their homogeneous analog. 

Hence, in the remaining part of this section, we turn back to the study of the homogeneous functional 
equation (14.22b . We shall in particular explain how one can infer the expected form of the Q-solutions to (14.22b 
from the quasi-periodicity properties of the coefficients, and we shall prove the completeness of these solutions 
for t(A) G S j in the case where N is even. 


4.3 On the reformulation of the SOV characterization of the spectrum in terms of solutions of 
an homogeneous functional T-Q equation: a few preliminary considerations 

Before turning to more precise considerations about the analytic properties of the Q-solutions to (14.22b . let us 
precisely formulate at which conditions the knowledge of some function Q( A) defines, through the equation 
(14.22b . a function t(A) which is an eigenvalue of the transfer matrix. 

6 Note that it is a priori possible to lower the degree M of the considered class of solutions Q{ A) by allowing some unknown 
parameters in /(A), as far as the total number of unknowns (i.e. of Bethe roots) is still at least equal to N. For instance, with the choice 
14.26b . one may consider an equation 14.25b for N unknowns which would be the set of the M — 1 roots of Q( A) and the parameter fi of 
/(A). We would then obtain Bethe equations of slightly different form coupling these two subsets of unknowns. 

7 In fact, if 77 £ R, they may be some additional restrictions on the set of parameters ft for which this assertion is valid (see 
Theorem IB. lb . Moreover, we also suppose that fj, is such that fi + (N — M)?? — £j, H + (N — M — 1)77 — £ 7 , fi — to,o — £.j, 
V + n- to ,o - & i r, V) € (1,..., N}. 
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Theorem 4.2. Let us suppose that the inhomogeneity parameters £ 1 ,..., £n °f the model satisfy (13.81) . and let 
t(A) be an entire function such that the following conditions are satisfied: 


1. there exists a function Q(X) such that t(A) can be expressed as 


z< \ \ — A x,y(^) Q(^ ~ v) + D (^) Q(^ + d) 

t(A) = m : 


(4.27) 


2. the function Q( A) of\7] is such that, for each n G {1,..., N}, (Q(£ n ), Qi^n ~ r l)) X (0,0); 

3. t(A) satisfies the quasi-periodicity properties (14.171) and (14.181) . 

Then t(A) is an eigenvalue of the K-twisted antiperiodic transfer matrix (i.e. t(A) G S^=), with corresponding 
eigenstates given in terms of the function Q( A) oj\ 7] as 


l^ K) ) 

= E 

n 


hG{0,l}^ 

1 0=1 



N 

<^ K) I 

- E 

n 


hG{0,lF 

1 0=1 


1 p iyn A x,y(Ca) \ n(c{ha)\ 

D (Za-V)J } 


det l"0(°' h )l |h,0), 
N L J 


(4.28) 


(4.29) 


As discussed in the previous subsection, this theorem is a direct corollary of the SOV characterization 
of the antiperiodic transfer matrix spectrum and eigenstates of Theorem 14.11 which can easily be proven by 
particularizing the relation (14.271) at the 2N points £ n , — q, 1 < n < N, so as to recover (14.211 ). The whole 

problem is therefore 

(a) to characterize the functional form of the function Q(X), preferably in a completely factorized form, so as 
to rewrite the entireness condition for t(A) as a system of Bethe equations for the roots of Q( A); 

(b) to prove that, for each t(A) G S T , there indeed exists such a Q-solution to the corresponding T-Q equa¬ 
tion, i.e. to prove the completeness of the aforementioned characterization in terms of solutions to Bethe 
equations. 

Before turning to these more delicate points that we shall partially solve in the next subsections, we would like 
to make a few remarks about the formulation of Theorem 14.21 which can in fact be rewritten in slightly different 
equivalent forms. For instance, it is easy to see that the condition 13 can be replaced by some equivalent 
conditions on the function Q(X): 

Proposition 4.1. For some given function Q( A), we define a function t(A) by the relation (14.271) . Then 
(i) the function t(A) defined by (14.271) satisfies the quasi-periodicity property (14.171) if and only if the function 
Wq\\) = Q(X + 7r) Q(X - n) - (-l)y Q(X + TT - n) Q(X) (4.30) 

satisfies the relation 

D (A) Wq\x + q) = (-l) x+xy a(A) Wq^ (A); (4.31) 
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(ii) the function t(A) defined by (14.271) satisfies the quasi-periodicity property (14.181) if and only if the function 


(A) = Q(X + ttuj) Q(X — 77 ) — (-l) x e"Q(A + ttw - 77 ) Q{ A) (4.32) 

satisfies the relation 

D(A)iy^ 2) (A + 7?) = (-l) y + xy e -^A(A) W^ 2 ) (A). (4.33) 

Hence the condition [?] in Theorem \4.2\ can be replaced by the condition: 

3’. the function Q(\) o/[7] satisfies the relations (14.311) and (14.331) . 

Another remark comes from the fact that, if Q(A) satisfies (14.221) for some function t(A) satisfying (14.171 ) 
and (14.181 ). then Q(A + it) and Q(X + 7 ru>) satisfy respectively the following equations: 

t(A) Q(A + 7r) = (-l) y A Xj y(A) Q( A + 7r — 77 ) + (-l) y d(A) Q( A + 77 + 77 ), 

t(A) Q( A + ttoj) = (-l) x e~ lNv A x>y (A) Q( A + 7 rw - 77 ) + (-l) x e* N,? d(A) Q(A + 7 ru; + 77 ), 

which means that we can in fact slightly relax the admissibility condition 12 ] according to the following propo¬ 
sition: 

Proposition 4.2. If Q( A) is a solution to (14.221 ) for some function t(A) satisfying (14.171 ) and (14.181) . then 
e ±,7ry 7 Q( A + 7 r), Q(A + 7ru;) am/ e*' N+ei 7“ +e2 7“^ A Q(A + 7r + ttlo) (with 61,62 = ±1) are also 

solutions of (14.221) for the same function t(A). Hence the condition \2\ in Theorem \4.2\ can be replaced by the 
condition: 


2’. the function Q( A) of\T\ is such that, for each n G {1,..., N}, there exists (a n ,/3 n ) G {0, l} 2 such that 
{Q{in + OLntt + fin™), Q{fi n + a n 7T + /3 n tTU - 7?)) (0, 0), 

with corresponding eigenstates given by 


i*i ) = e n 

hG{0,l} N o=l - 


( ei K ?+ «N + =,id*-> Q(e .) +ttoir + ^ u) 

V KD (Co -mJ 


x det [0(°’ h) ] |h,0), (4.34) 


(47^1= [] + det[0(°’ h )](O,h|, 

he{o,i} N o=t L J 

(4.35) 


t(^a) 


instead of (14.281) . (14.291) . 


Note finally that the relation (14.311) (respectively the relation (14.331) 1 can be understood as some “wronskian- 
type” identity for the two solutions Q( A) and e ±,7ry 7 Q (A + n) (respectively Q{ A) and e ^ N± 7“^ A Q (A + 7rw)) 
of the functional T-Q equation (14.221) . as a particular case of the following very general property: 
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Proposition 4.3. Let Qi(A) and Q 2 (A) be two solutions of the functional equation (14.221) for some given function 
t(A). Then their quantum Wronskian, 

W 12 (\) = Qi(A - rj) Q 2 (A) - Qi(A) Q 2 (A - 77 ), (4.36) 

satisfies the relation 

d(A) Wi 2 (A + 77 ) = (-l) x+y+xy A(A) Wi 2 (A). (4.37) 

Hence 

W 12 (X) = / 12 (A) d(A), with / 12 (A + 77 ) = (-l) x+y+xy / 12 (A). (4.38) 

Remark 4.1. If Q(A) satisfies the homogeneous functional equation (14.221) then, for each fc € Z, e 2 ” fc 7 Q(A) 
also provides a solution to the functional (14.221 ). which is however not independent from Q( A) in the sense that 
then quantum Wronskian (14.361) is identically zero. 


4.4 Study of the homogeneous T-Q functional equation: an Ansatz for the Q-solutions 


As iust announced, the neriodicitv nronerties of the coefficients of (14.221) for t(A) G s r enable one to make a 
reasonable guess about the functional form of its possible entire Q-solutions, i.e. the solutions that are suscep¬ 
tible to lead to a system of Bethe equations. The idea is to transform this equation, similarly as in E71 . into a 
difference equation with elliptic coefficients: 


Q(A + -,) + /, (A) 7 222L Q( A) + A(A) LLftl Q( A - ,) = 0. 
Here we have set 


hi A) 


t(A) 
d (A) ’ 


/2(A) 


Ax,y(A) 

D(A) 


and Q( A) 


Qi A) 

p(\y 


(4.39) 


(4.40) 


where P( A) has to be chosen appropriately so as to ensure that the coefficients of (14.391 ) are doubly-periodic. 
Since the functions /,(A) exhibit different quasi-periodicity properties according to the values of x and y, we 
shall now consider the three different cases x = 0, y = 0 and x = y separately. To this aim, it is convenient to 
introduce shorthand notations for different variants of the theta function that we shall use within the study of 
these different cases. We therefore define the following odd functions of A, 


<wa)^4!|) 


# y=0 (A) = $i(A | 2ui), 
X 
.2 


6 x=y (X) = eh 6 t 


u) 01 


A + 7T + ITU) 


U) 


which satisfy respectively the quasi-periodicity properties: 


(4.41) 

(4.42) 

(4.43) 


2 

0 

II 

X 

1 

0 x =o(A + ttuj) = —e lX 0 x =o(A), 

(4.44) 

2 

0 

II 

0 y =o(A + 2itlo) = _ e - 2iA - 2i ™ 0 y=o (A), 

(4.45) 

2 

X 

II 

X 

1 

0 x=y (A + 2t rw) = - e ~ 2iX ~ 2i ™ 0 x=y (A), 

(4.46) 

-) = ±ie- lA -^0 x=y (A). 


(4.47) 
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(4.48) 


The case x = 0. Let us choose P( A) in the form 

N 

P(A) = []^x=o(A-^), 

3= 1 

for arbitrary Zj, 1 < j < N. Then the coefficients /i(A) and / 2 (A) p|^ + ^| for the equation (14.391) for 

Q( A) are elliptic functions of periods 2ir and 7 tui. Following li27ll . we may therefore look for the double-Bloch 
solutions of (14.39b . i.e. for the meromorphic solutions Q( A) = such that 

Q{\ + 2n) = B\ Q(A), and Q{\ + too) = B 2 Q(\), (4.49) 


for some Bloch multipliers B\, B 2 . It is easy to see l28l that any double-Bloch function of the form (14.491 ) can 
be written in the form 


Q(A) 


M 

c 0 e" A n 


3 =1 


^X=0(A ~ A j) 

#X = o(A /ij) 


(4.50) 


for some integer M, some complex parameter a and some sets of zeroes and poles A j and fij, I < j < M. If 
Q(A) of the form (14.501 ) is solution of (14.391) . then Q( A) = P( A) Q( A) is solution of (14.221) . and we expect the 
latter to be cntirtQ. which means that it can be written in the form 
N 

Q(\) = c Q e QX l{0 x= o(\-\j), (4.51) 

3= 1 


for some complex parameter a and some set of roots Xj, 1 < j < N. 

Note that, in the obtention of the Ansatz (14.511) for Q(A), we have only used a weaker periodicity property 
than (14.171) for t(A). Imposing that t(A) strictly satisfies (14.171) is equivalent to the additional constraint (14.311 ) 
for the corresponding function Wq(X) defined in term of Q(A) as in (14.301 ). It means that 

Wq 1} (A) = 5 (1) (A)d(A), (4.52) 

where g W (A) is 77-periodic: (A + ij) = g^ (A). On the other hand, for Q(A) of the form (14.511) . it is easy to 

see that Wq^ (A) satisfies the quasi-periodicity properties 

Wq\\ + tt) = (-l) N +y +1 e 2 ™ ^(A), (4.53) 

w §\A + 7 TU) = (_ e -2^-^) N e 2iJ2X j +M V +2nua W £\\), (4.54) 

which have to be compared to the quasi-periodicity properties of d(A). Hence g ll >(\) also satisfies the quasi¬ 
periodicity properties: 


S f ( 1 ) (A + 7r) = (—l) y+1 e 2na (4.55) 

3 (1) ( A + ttw) = e 2< [i: A J+5’»-E^]+ 2,rw “ & (i)(A). (4.56) 


8 It is easy to see that, if Q( A) solution to 14.221 admits some pole g,j, then it should have an infinite number of poles of the 
form /ij + kg for an infinite number of k £ Z, which is clearly not compatible with the form 14.501 of Q( A) (we suppose here that 
g (/z 7tQ + nujQ, and that the inhomogeneity parameters satisfy 13.81 ). 

9 Indeed, if Q( A) of the form 14.511 is a solution to 14.221 . e l ~ x Q( A + 7r) is another solution to 14.221 that we expect here to be 
independent from the first one. 
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(4.57) 


It means that there exist two constants c^J and a i such that 

0 ( 1 ) (A) = c$e“ lA , 


and that there exist k \, k 2 ■ k :i G Z such that 
a\r] = 2niki } 

a\i\ = *7r(l — y) + 27ra + 2irik 2 , 


aiiruj = 2 i 




+ 27ru;a + 2-Kik^. 


Hence, if c^J / 0, one obtains the following sum rules: 


a = i 

N 


y - 1 , , * ' 

—^- Vk i- k 2 

2 T] 


^2 A i = - IT 7 ? + (! - y)-7T + fe 2 vru; - fc 3 7 r, 


, j=i fc=i 

for some ki,k 2 ,k 3 G Z. 


(4.58) 


(4.59) 


Note that, if there indeed exists such a solution, it means that there exists such a solution for which k 2 = 0 
and k\, k$ G {0,1} (this is due to Remark [4711 and to the quasi-periodicity properties of the function 0 x =o with 
respect to a shift of one the roots A j by 27r or ttcj). In fact, one can even be more precise, and formulate the 
following Ansatz: 


• if (x, y) = (0,1), we expect two independent solutions of the form 

N 

Q(A) = n 0 x=o(A - Xj), 

3= 1 

with 

N N N 

= ^2^ k -—’l + k-K, k G {0,1}, 
j =1 k =1 

and 

Q(A) = e^ A Q(A + vr); 

• if (x, y) = (0. 0) and N odd, we expect two independent solutions of the form 

Ik 1 N 

Q(A) = e HT-iJ A JJ 0 X=O (A — Aj), k G {0,1}, 

3 = 1 

with 

N N M 

E v —' N 7 TLO 

x 3 = 2^^-2 71 + ~T’ 

3 =1 fc=l 

and 

Q( A) = Q(A + 7r). 


(4.60) 


(4-61) 


(4.62) 


(4.63) 


(4.64) 


(4.65) 
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(4.66) 


The case y = 0. A similar reasoning can be made in this case, by choosing P( A) of the form 
N 

P(A) = n^y=o(A-^), 

7=1 

with arbitrary Zj, 1 < j < N, so as to obtain elliptic coefficients with periods ir and 2ttco in (14.391) . It leads to 
the following Ansatz for Q(A): 

N 

Q(A) = cq e aX ]^[ 0 y= o(A — Xj), (4.67) 

3 = 1 

for some complex parameter a and some set of roots \j, 1 < j < N. The condition (14.331) . which is a necessary 
and sufficient condition for the corresponding function (14.27b to satisfy (14.181) . results into the following sum 
rules for (14.671 ): 


a = i 

N 


h--k 2 

v 

N 


^ Aj = - ^ + (i - x)| + k 2 iru - k 3 7T, 

„ 3 =1 fc=l 


(4.68) 


for some ki,k 2 ,k 3 G Z. 


As previously, we can be more precise. Indeed, if such a solution exists, then it is easy to see, from Re- 
mark l4.1l and considerations about the quasi-periodicity properties of the function 9 y= o with respect to a shift of 
one of the roots \j by 7r or 2Tree, that there also exists a solution for which k 3 = 0 and &q, k 2 G {0,1}. Then 
one can formulate the following Ansatz: 


if (x, y) = (1,0), we expect two independent solutions of the form 
N 

Q(A) = e~ ikX n 9 y= o(A - Xj), k € {0,1}, 

3 = 1 


with 


N N 


N 


J2 x i = ~ v 1+ /c7rw ’ 


3 = 1 k=l 

and 

Q(A) = e * (N+ 7^ A Q (A + 7 rw); 

if (x, y) = (0, 0) and N odd, we expect two independent solutions of the form 
N 

Q(A) = e ii T A JJ 6» y=0 (A - Aj), A: € {0,1}, 

7=1 


with 


N N 

E \ ^ N 7T 


7=1 fe=i 


and 


Q(A) = e* NA Q(A + 7ra;). 


(4.69) 

(4.70) 

(4.71) 

(4.72) 

(4.73) 

(4.74) 
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The case x = y. In that case, the choice 


N 

PW = ]\^ =y (\-z 3 ), (4.75) 

3 = 1 

with arbitrary Zj, 1 < j < N, leads to an equation (14.39b which has elliptic coefficients with periods itlo + i r 
and 7 tui — ir, and to the following Ansatz for the solution Q(X) of (14.221) : 


N 

Q(A) = c Q e“ A n^x= y (A-A i ), 

3 = 1 


(4.76) 


for some complex parameter a and some set of roots Xj, 1 < j < N. It is then easy to see that the conditions 
(14.311) and (14.331) result into the saniJ^I sum rules for the parameters entering (14.761) . which can be written as 


a = i 

N 


y - 1 , , Ti- ' 

—--b «1- k-2 

2 rj 


k =t 


J2 A i = + (tt + 

L 3 = 1 


7 TUJ 


i — y 


+ k 2 ) + k 3 n, 


(4.77) 


for k \, k 2 , k 3 G Z. 

Here again one can be more precise: if such a solution exists, it is easy to see that, due to Remark |4~TI and 
considerations about the quasi-periodicity properties of the function 9 X=V with respect to a shift of one of the 
roots A j by 7T + Tree or 27r, it exists for k 2 = 0 and k \. k 3 G {0,1}. Due to the existence of a second solution of 
the form 


i —A \ | \ 

e * Q {A + 7t) oc e 


*(N+f)A 


Q(A + 7rw) 


(4.78) 


one can finally formulate the following Ansatz: 


• if (x, y) = (1,1), we expect two independent solutions of the form 
N 

Q(A) = nM^-Aj). (4-79) 

3 = 1 

with 

N N N 

Aj = + fe7T, k G {0,1}, (4.80) 

7=1 fc=i 

and 

Q(A) = e*7 A Q(A + 7r) oc e 4 N+ 7) A Q(A + 7ro;); (4.81) 

10 In fact, for Q(A) of the form <476t . one has tk^ 2) (A + n) = (-l) x+1 (ie~ iX - il! ?) N e iE *i+A+™0« W«(A). 
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if (x, y) = (0, 0) and N odd, we expect two independent solutions of the form 


Ik il N 

Q(A) =e i l ! f- 1 ii X k€ {0,1}, 

3 =1 

with 

N 7T + 7 TOJ 

j=1 k=1 

and 


(4.82) 


(4.83) 


Q(A) = Q(A + ^) oc e ,:NA Q(A + 7ro;). (4.84) 

Remark 4.2. The case (x, y) = (0, 0) with N odd can be considered from three different viewpoints (as a 
particular case of x = 0, of y = 0 or of x = y), hence leading to three possible different systems of independent 
solutions ( (14.631) and (14.651) . (14.72b and (14.74b . or (14.82b and (14.84b ). 

In the previous study, we have exhibited possible forms for the entire Q-solutions of (14.22b . hence suscep¬ 
tible to lead to a system of Bethe-type equations. One can in fact be slightly more precise in our argument by 
remarking that, if Q{\) is a solution of the form (14.51b . (14.67b or (14.76b to the homogeneous functional equation 
(14.22b for some function t(A) satisfying (14.17b and (14.18b . then at least one of the two functions (14.30b or (14.32b 
should be non identically vanishing. This means on the one hand that the sum rules (14.59b . (14.68b or (14.77b 
should indeed be satisfied, so that we can indeed refine the form of Q(X) as in (14.60b . (14.63b . (14.69b . ( 14.72b . 
(14.79b or (14.82b according to the case. This means on the other hand that the condition 2’. of Proposition I4.2l is 
automatically satisfied. This is due to the following proposition and corollary: 

Proposition 4.4. lfQ{ A) is a non-zero solution of (14.22b of the form (14.51b . (14.67b or (14.76b to (14.22b for some 
function t(A) satisfying (14.17b and (14.18b . then at least one of the two functions Wq\ A) (14.30b and Wq\ A) 
(14.32b is not identically zero. 


Proof Let us suppose that both functions Wq - (X) and Wq\ A) are identically zero, i.e. that Q(X) satisfies the 
relations 




Q{ A + 7r) 

Q{ A + 7r — rj) 
Q{ A + 7 tuj) 
Q (A + 7 tu; - rj) 


(-1 ) y 


Q{ A) 

Q {A - n) ’ 


(_l) x e-Mr, 


Q( A) 

Q(X-r)Y 


(4.85) 


This means that the function 


Q{X-ri) ' s a double-Bloch function which can therefore be written in the form 


Q(A) 
Q (A - rj) 


cp e 


/3A 


M 

n 

3 = 1 


g (A ~ Pj) 

6 (A - uj) 


(4.86) 


for some integer M, some complex parameter j3 and some sets of zeroes and poles fjj and vj, 1 < j < M. 
It follows that Q( A) should be of the form (14.23b which, as already mentioned, is not compatible with the 
quasi-periodicity properties of (14.22b . □ 
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Corollary 4.1. Let t(A) satisfying (14. 171) and (14. 181) and let us suppose that there exists a function Q(X) solution 
to the homogeneous functional equation (14.221) associated with t(A). 

• If x = 0 and if Q{ A) is of the form (14.511) . then Wq\ A) is not identically zero and the sum rule (14.591 ) is 

satisfied. Moreover, for all n £ {1,... N}, (Q(Cn — Q(£n — ?? + vr)) ^ ( 0 , 0 ). 

/r)\ 

• If y = 0 and if Q{X) is of the form (14.671) . then Wq (A) is not identically zero and the sum rule (14.681) is 

satisfied. Moreover, for all n £ {1,... N}, (Q(£, n ~ v)i Q(^n ~ V + TTUj)) / ( 0 , 0 ). 

• Ifx = y = 0 and if Q(X) is of the form (14.761) . then Wq\\) and Wq\ A) are not identically zero and the 
sum rule (14.771) is satisfied. Moreover, for all n £ {1,... N}, ( Q(^ n ~ ??), Q(£n — p + 7r)) f (0,0) and 

{Q{in ~ ?7), Q(£n - r] + ttcj)) / (0,0). 

Moreover, in all these cases, t(A) is a transfer matrix eigenvalue. 

Proof. Let us first consider the case x = 0. The first part of the assertion is a direct consequence of Propo¬ 
sition |44j of the fact that Wq\ A) is identically zero for Q(X) of the form (14.511) and of the previous study. 
Therefore Wq ' 1 (X) = cjy e° |A t)(A) with c^- f 0. This is clearly incompatible with the possibility that, for 
some n £ {1,..., N}, ( Q(£n — d)-> Q{Cn ~ ?? + 7 r)) = (0, 0), since in that case we would have Wq (£n~ rj) = 0. 
Hence the second assertion follows. 

The cases y = 0 and x = y can be proven similarly. □ 

It is also easy to see that, whenever such solutions exist for t(A) £ they are unique, which means that 
the characterization of a given eigenvalue (and corresponding eigenvector) of the antiperiodic transfer matrix in 
terms of the corresponding type of solutions of Bethe equation is uniquely determined (i.e the same eigenvalue 
cannot be obtained by two different solutions of the same type). 

Proposition 4.5. For each t(A) £ the equation (14.221) admits 

• at most one independent solution of the form (14.601) if(x, y) = (0,1); 

• at most one independent solution of the form (14.691) if (x. y) = ( 1 ,0); 

• at most one independent solution of the form (14.791 ) if (x, y) = (1,1); 

• at most one independent solution of the form (14.631) with the constraint (14.64b . at most one independent 
solution of the form (14.721) . and at most one independent solution of the form (14.821) with the constraint 

<ra>r(x,y) = (o,o). 

Proof. Let us suppose that there exists, for (x, y) f (0,0) and for a given t(A) € two solutions (f (A) 
and CfiiX) of the form (14.601 ) (respectively (14.69b . respectively (14.79b ). Then their quantum Wronskian (14.361 ) 
satisfies the periodicity property 1 L| 2 (A + 27r) = lL"i 2 (A). Hence !L| 2 (A) = / 12 (A) d(A), where / 12 (A) is an 
entire function with two incommensurable periods 27r and 2 // (we suppose that // / vrQ), i.e. it is a constant. 
The latter is equal to zero if (x, y) 7 ! (0, 0) due to (14.381) . 

Let now (x, y) = (0, 0) and let us suppose that, for a given t(A) € there exist two solutions Q 1 (A) and 
Q 2 (A) of the form (14.631 ) with the constraint (14.641 ) (respectively of the form (14.721) . respectively of the form 
(14.821 ) with the constraint (14.8310 . By considering the quasi-periodicity properties of their quantum Wronskian 
(14.361) 1 L| 2 (A) = / 12 (A) d(A) with respect to shifts of 2tt and iru (respectively of it and 2 Tree, respectively of 27 t 
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and ±7r + 7 tu), one obtains that / 12 (A) is of the form / 12 (A) = ci 2 e“ 12A , and one obtains some conditions on 
a 12 if C 12 f 0. These conditions are clearly incompatible with the supposed form of Qi(A) and Q 2 (A), which 
means here also that C 12 = 0 and that the two solutions are not independent. □ 

We shall see in the next subsection that the reformulation of the spectral problem for T(A) in terms of 
solutions of the form (14.601) . (14.691) or (14.791) is effectively equivalent to the SOV characterization of the spectrum 
of Theorem 14. 11 at least in the case N even. 

4.5 On the reformulation of the SOV characterization of the spectrum in terms of solutions of 
an homogeneous functional T-Q equation: proof of the completeness for N even 

We shall now prove the converse of Theorem 14. 2 1 in the case of even N and for Q( A) of the form (14.601 ). (14.691 ) 
or (14.791) . i.e. the completeness of the solutions to the corresponding Bethe equations. More precisely, we shall 
prove the following result: 

Theorem 4.3. Let us suppose that the inhomogeneity parameters £ 1 ,... ,£n satisfy (13.81) and let t(A) be an 
eigenvalue of the antiperiodic transfer matrix T(A) (i.e. t(A) € Xy). Then, if N is even, there exists a function 
Q{ A) of the form 

N 

Q(A) = Y[0 x (\-\ j ), (4.87) 

5 = 1 

for some set of roots Ai,..., A|\|, such that t(A) and Q( A) satisfy the homogeneous functional equation (14.221) . 
Here the notation #x(A) stands, respectively, for the function (14.411) if (x, y) = (0,1), for the function (14.421) if 
(x,y) = (1,0), or for the function (14.431) if (x, y) = (1,1). 

From Proposition 14.51 we know that such a solution is unique (up to normalization). We also know from 
Corollary 14.11 that this solution is such that condition 2 ’. of Proposition 14.21 is satisfied. Hence, in the case 
of even N, it means that the discrete SOV characterization of the spectrum is equivalent to the description in 
terms of Bethe equations based on solutions of the form (14.871) of the homogeneous functional equation (14.221) . 
The description of the eigenstates can then be obtained as in Theorem 14.21 or more generally as in its variant 
Propostion 14.21 Note here that it is possible to rewrite the corresponding eigenvectors in a form more similar 
to what we have in the ABA framework, i.e. by multiple action, on a given reference state, of a product of 
operators evaluated at the corresponding Bethe roots. However, due to the fact that the function Q{ A) (14.871 ) 
is not a theta function of the same type (i.e. with the same quasi-periods) as the other functions defining the 
model, and unlike what happens in the context of the reformulation of Appendix |Bj (see Corollarv lB.il) . we do 
not use the whole operator V to construct the eigenstates, but only “parts” of this operator: in the present case, 
the operator V is split into a product of commuting operators, and these are the latter which are used to construct 
the eigenstates. More precisely, for each N-tuple (3 = (/3i,..., /3 n) £ {0,1} N , let us define in End(D(6VD),i\i) 
the following diagonal operator on the SOV basis given by the vectors |h, r ): 

N 

h,r) = Y\ 

71=1 



*7r(x+y— xy)h n 

l N 


+is y= 0 (\-d hri) ) 


#x(A - £i hn) + /3 n vr x ) l |h, r), 


(4.88) 
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where 


c x — < 


04(O|W) 

ie^TT 0 2 (O|w 

1 
2 

7T X = (1 - 5 y=0 ) 7T + (5 y=0 7TW. 
We have that 


he *2 0 2 (O|w) 0 3 (O|w) 0 4 (O|w) 


if x = 0, 
if y = o, 
if x = y, 


t7r(x+y—xy) \-^N 


(4.89) 

(4.90) 

(4.91) 

(4.92) 

(4.93) 


Dp{X) D 1 _ /3 (X) |h, r) = e —JJ 0(A-^) \h,r), 

71 = 1 

where 1 — /3 stands for the N-tuple (1 — /3i,..., 1 — /3 n)- This means that 

P(A) = T d £> /3 (A)£>i_ /3 (A), 

where To is the shift operator defined on the SOV basis as 

_ . t7r(x+y-xy) s h -s 0 • Sh“ s l 0(t,-+1Vi) 

T d | h, r) = e N 2 2 -A-Lth-L |h, r + 1). 

"(tr+1,1) 

Then one can formulate the following corollary: 

Corollary 4.2. Under the condition (13.81 ) and if N is even, there exists a one-to-one correspondence between 
Tiy and the set T^bae of different (up to the real quasi-period of the function Ox) Bethe roots A = {Ai,..., An} 
such that 

1. the function 

. Q( x ~v) , Q( x + d) 

A ^ (A)_ g(Ar + ()_ w { 1 

is entire and satisfies the quasi-periodicity properties ( 14. 1 7I )-( 14. 1 81 ). 

2. Vn € {1,..., N}, there exists fi n E {0,1} such that Q(£ n + /3 n Tr x ) f 0, 

where Of A) is defined in terms of A bv (14.871). The eigenvalue t(A) ee t associated with A E 'Fhae is then 
given by the entire function (14.941 ). For any k E C \ {0}, the corresponding one-dimensional right and left 
eigenspaces of the n-twisted transfer matrix T^\ A) are the one-dimensional subspaces spanned 

by all vectors of the type 

N N 

lO = II D ^ X i) respectively (T^| = (f2^| D p {\j), (4.95) 

i = 1 i =1 

for any N-tuple (3 E {0, 1 } N . In (14.951) . the operators Dp (A) are defined as in (14.881 ). and the reference states 
\n^) and (n(*)| are 

l« (K) >= £ n(S(|^)'‘“^[e<»' h »]|h,o>, ,4.96) 


he{0,l} N 0=1 

N h 

(o( k )| = n ( Ke * y,? ) a det [ 0(o,h) ] (°> h i- 

he{o,t} N «=i 


(4.97) 
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The proof of Theorem 14.31 is based on the following lemma: 

Lemma 4.1. Let us suppose that the inhomogeneity parameters £ 1 ,... , £n satisfy (13.81) and let N = 2M be 
even. Then, for each t(A) £ YLj, there exist A £ C and two non-zero entire functions Q+( A) and Q-( A) with 
the following quasi-periodicity properties: 


Q ± (A + vr) = (±l) y Q ± (A), 

Q±(A + nut) = (±l) x ( — e —: 

which satisfy the following system of 2N equations: 

f t(0) Q+(Zj) = -a(£, 0 Q-(0 - ??), i < . 

1 t(£? - ??) ~v) = d(0 - ??) Q+(&), 


(4.98) 

(4.99) 


(4.100) 


Proof On the one hand, since t(A) £ the system (14.1001 ) is equivalent to the following system of only N 
equations: 

t&) Q+&) = -A(Ci) Q-& -y), i <j< N. (4.101) 

On the other hand, the entireness and quasi-periodicity properties of the two functions Q±(A) are equivalent to 
the fact that the functions 

Q±( A) = e Q ± ( A) (4.102) 

arc theta functions of A of quasi-periods ( 7 r, 7 rut), of order M and of norm 


A± = A- 


1 =F 1 / 7r ttuj , 

V y ir>- 


This is also equivalent to the fact that the functions Q±(A) can be represented in the following form. 


M 


0+<A) - E e «c£ 1 ft,-A +) n f-y «♦«..>, 


(4.103) 


(4.104) 


fc=i 

M 


l^k 


Q- (A) = £ 


k=1 


,i(y-M)(A-€ ifc +?j) ^ ~ + y + ~ Mp - A_) 

xn ^ r «-&-* (4105 > 


1=0 

e^k 


in terms of some arbitrary M-tuples (<3 + (£q),..., <3+(£i M )) and (Q-(& - 77 ),... , Q-(£i M - 77 )), provided 
that A + — ^ I and A_ — + M// f T. Here, we have arbitrarily split the set of the N 

inhomogeneity parameters {£ 1 ,... ,£n} into two disjoint subsets {£j 1; ...,£,; M } and {£j M+1 ,... ,£« N } of the 
same cardinality M = N/2. 

Hence, the system (14.1011 ) for £)+(A) and Q-( A) is equivalent to the following linear system for these two 
M-tuples: 

«-(&,- 7 ) = «+(&,), 

M j = 1, . . . , M, 

E[^(A )] jk Q+& k ) = 0, 

, k= 1 


(4.106) 
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where X\ (A) is the M x M matrix of elements: 


- TT _4) 


[^(A)] afe = e- lM ^ + M-^)n 


/ = 1 0 tei“&J 

l^b 


^tea + M £*& + £^=lte* A+) t(^*a + M) 


0 (£^-Ah 


Ate, 


a+M7 


— e *y(Si a+M 


-Zi b ) ^tea + M C» b + El=l &t ^ 7 / A-) ttefc) 

^(Efcl^-Mry-A.) AfeJ 


(4.107) 


The first line in (14.1061 ) corresponds to the equations i\,... ,im of (14.1011) . which completely fixe the M-tuple 
(Q-tei - T]), ■ ■ ■, Q-te M “ V)) in terms of the M-tuple (Q+tei), • ■ ■, Q+te M ))> whereas the second line 
in (14. 1 06b corresponds to the equations *m+i, ■ ■ • ,*n °f (14.1011) . in which we have used the representations 
(14. 1041) and (14. 1051) that we have rewritten (by means of the first line of (14.1061 0 in terms of the M-tuple 

( Q + tei )) ■ • ■ . Q + teiJ ) onl y - 

Hence, this system admits a non-zero solution if and only if the determinant of the matrix A*(A) (14.1071 ) 
vanishes. It is easy to see that this determinant is a quasi-periodic function of A with quasi-periodicity properties 


det [AHA + 7r)l = det [AHA)]. 

h/l L J M L J 


(4.108) 

M L '' ' J M L " ' J ' 

det [AHA + 7rw)l = e 2l ^ a =^ hqc(A)], (4.109) 

M M L J 

so that, if not identically zero, it is a non-co ns tan 0 function of A which can be written in the form 

n^A-A,) 


det [At(A)] = cx - M M 

M L J 0(A + -£, M = 1 ^) M 0(A_-£, M = i e^ + Mr / ) M 


(4.110) 


in terms of some roots A j which are not all equal to the roots appearing in the denominator. This ends the proof 


of Lemma |4~T1 


□ 


Proof of Theorem 14. .4 Let t(A) € S T and N = 2M be even. Then Lemma |4~TI implies that there exists 
two non-zero entire functions Q±(A) which satisfy the quasi-periodicity properties (I4.98I) - (I4.99I) and the system 
(14.1001 ) for some A G C. This implies that there exist two entire functions f±( A), with quasi-periodicity 
properties 


/ ± (A + vr) = (±l) y /±(A), (4.111) 

/±(A + 7tw) = {±l) x (- e - 2iX ~ in0J ) M e 2iA ~ m ™ f±{\), (4.112) 


which satisfy together with the functions Q±(A) the following system of functional equations: 

r t(A) Q+(A) = —a(A) Q_(A - q) + d(A) /_(A + v ), 

\ t(A) Q-(A) = — a(A) /+(A — rj) + d(A) Q+(X + rj). 

Particularizing the first line of (14. 1131) at the points £j — rj, j = 1,..., N, and the second line at the points £j, 
j = 1,.... N, and using the fact that t(A) £ we obtain that 

<2+te ~v)Q-(€j) = /+(0 1 <J<N. (4.114) 

n at least for some adequate splitting of (£i,... , £n}- 
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Note that the function Q + (X — //) Q-{ A) and the function f+(\ — >]) /-(A) satisfy the same quasi-periodicity 
properties: 


F-\ —(A + 7 r) = -F-i—(A), (4.115) 

F_|_(A + 7TU;) = (_ e - 2iA - i,raj ) N e 4?;A+iN '?- lN ^ F+_(A), (4.116) 

so that F h _(A) is a theta function of order N and of norm 2A + Mp — Mtto;, where F. ( _(A) stands either for 

the function Q+(A — 77 ) Q_(A) or for the function /+(A — rj) /_(A). Hence (14.114b implies that the identity is 
in fact valid at the functional level, i.e. for any A € C: 

Q + (A - r?) Q_(A) = /+(A - 77 ) /-(A). (4.117) 

One can therefore eliminate / + (A — 77 ) in the system (14.113b . and we obtain the following functional system for 
the function /_(A): 

' t(A) Q+( A) = -a(A) Q_(A - rj) + d(A) /_(A + p), 

( t(A) Q_( A) = — a(A) Q+(A ~^ ) Q ~ (A) + d(A) Q + (A + 77 ). (4 ' 118) 

Note that the ratio in the second line of (14.118b is in fact an entire function due to (14.117b . This systems implies 
that 


a(A)/ 1 (A) = d(A)/ 2 (A + t ? ), 


(4.119) 


in which we have set 

/t( A) = Q-( A - v) Q-( A) - Q+(A ~(I) Q " (A) Q+(A) ’ (4 ‘ 120) 

/ 2 (A) = /-(A) Q_(A - r?) - Q+(A) Q+(A - r?). (4.121) 

/i(A) and / 2 (A) are two entire functions of A which are both theta functions of order N. and therefore the 
relation (14.119b implies that there exists two constants c\ , c 2 G C such that 

'/t(A) = Cl D(A), 

( / 2 (A + 77 ) = c 2 a(A), i.e. /i(A) = / 2 (A) = ci di 

Hence 

/-(A) = Q-(X), and /+(A) = Q+(A), 

which means that Q_(A) and Q+(A) satisfy the following system 

t(A) Q+(A) = —a(A) Q_(A - r,) + d(A) Q_(A + 77), 
t(A) Q_(A) = —a(A) Q + { A - 77 ) + d(A) Q + (A + 77 ). 


(A). (4.122) 


(4.123) 


(4.124) 


Let us now define the entire functions 

m= Q+W+ t Q - W . -d q (a)= Q+W-Q-W 


(4.125) 


34 


















with quasi-periodicity properties 


(4.126) 

(4.127) 


Q(X + it) — (i — y) Q(X) + y Q(X), 

Q{ A + 7 tuj) = (-e- 2a -^)M e 2iA ~ iM ™ [(1 - x) Q( A) + xQ(A)] , 

and 

Q(A + tt) = (1 - y) Q(A) + y Q(A), (4.128) 

Q(A + tto;) = (- e -»-“)Me 2iA -‘M™ [(1 -x)Q(A) +xQ(A)] . (4.129) 

From these quasi-periodicity relations, it is easy to see that Q(A) and Q( A) are both (up to a constant nor¬ 
malization factor) of the form (14.871) . Moreover they satisfy the following respective homogeneous functional 
equations: 

t(A) Q(X) = —a(A) Q{ A — t]) + d(A) Q(X + rj), (4.130) 

t(A) Q(X) = a(A) Q(X — rj)~ d(A) Q(X + rj). (4.131) 

This ends the proof of Theorem 14.3 1 □ 


5 Local operators and dynamical Yang-Baxter algebra 


In the next section, we shall compute determinant representations for form factors of local operators in the 
1 (6VD) n" basis of the eigenstates of the re-twisted transfer matrix. As in the algebraic Bethe Ansatz framework 
[26, 30], such representations are based on the solution of the quantum inverse problem, i.e. on the fact that 
one can reconstruct the local operators we consider in terms of the generators of the Yang-Baxter algebra. A 
particularly crucial point in this respect comes from the fact that the positions on the lattice of these local 
operators are given in terms of propagators written as products of transfer matrices, so that their action on the 
corresponding eigenstates merely contributes as simple numerical coefficients. 

For the dynamical 6-vertex model that we consider in this paper, the local operators of interest are essentially 
of two types: local spin operators and local height operators. In our framework, local spin operators correspond 
to elementary operators acting non-trivially on only one factor of the space tensor product Vn = , V n and 

can be expressed in terms of the usual basis of elementary matrices {En ; i,j <E {+,—}, 1 < n < N}of 
End(V|\i) defined by 


El 3 


( S +X 6 +,s < 4 ,* 5 -, A 

V~X S +’j < 5-4 < 5 - 7 / 


i,j G {+,-},! < n < N, 


(5.1) 


whereas the local height operators are dynamical operators acting non-trivially on the dynamical space D. In 
the physical context of classical SOS face models, the local operators that we should especially consider for the 
computation of correlation functions are the following (see for instance Section 3 of 11321 1: the spin operators 
E+ + and E~~, which are associated with the values of the classical spin variable (+ or —) on the n-th vertical 
bond of a given vertical line of the lattice, and the local height operators 6^\ t € f 0 + r/Z, which are associated 
with the values t of the local height variable on the j-th site of a given vertical line of the lattice. Note in 
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particular that the local height operators 6) ' at site 1 (the reference site) have the following simple action on the 
dynamical-spin basis (12.311) : 

( < 8>n=l( n > h n \) <g> (f(a)| = 5 t ,t{a) (®n=l( n ,h n \) ® (t(a)|, (5.2) 

(®n=ll n X)) ® l*( a )) = S t,t(a) (®n= i\n,K)) ® |*(a)). (5.3) 

Its action on the SOV basis (I3.4I) - (I3.5I) hence simply follows from (13.61 ): 


(r,h 


? 1} = 6, 


't t 


, h ( r M 


xt 1 ) 


h. 


= 4 


Ij. j. 


|h,r). 


(5.4) 


In the following, we shall denote this operator as 6) ; = c)/(r). 

Reconstruction of local spin operators in terms of the elements of the dynamical periodic monodromy matrix 
(12.111) has been obtained in l30l . generalizing to the dynamical case the proof of 1261 [Ml . In this paper, we need 
instead to express these operators in terms of the elements of the K-twisted antiperiodic monodromy matrix 
(12.201) (using in particular propagators given in terms of the rc-twisted transfer matrix (12.251) 1 so as to be able to 
easily compute their action on the eigenstates (14.3I) - (14.4I) . 

So as to remain as general as possible and to present a proof that is valid in both periodic and antiperiodic 
cases, we shall instead consider a quasi-periodic monodromy matrix of the form 


M ( 0 r) (A) = Y 0 M 0 (A), 


(5.5) 


where A4 is given by (12.1 II) and Y stands for any numerical invertible matrix (the ^-twisted antiperiodic case 
corresponding to the particular choice Y = X (f;, a x ). The corresponding quasi-periodic transfer matrix will be 
denoted by T (i )(A) = tro ^(A)]. The aim is therefore to express local operators in terms of the operator 
entries of (15.51 ). As in the periodic case, this reconstruction is based on the following lemma: 


Lemma 5.1. 001/ We have the following identity between products of monodromy matrices: 


k =l 


k -1 


nM afci i..jvKfc|T + r ? ^ 


Z=1 


% I = 


k -1 


a 2 ...afc_i ( £fc| r + V E "f+ i E 

k =1 '■ 1=1 l=k +1 


a. 


<H / ’ 


(5.6) 


in which the symbols (respectively ) means that the product is ordered from 1 to m (respectively from m 

k =1 k =1 

to 1). In this expression, M a fcfc+i..jv’ai a 2 ...a fc _ x denotes the monodromy matrix of a chain ofN sites labelled in 
this order by k , k + 1,... N, a\, 0 , 2 , ■ ■ ■ dk-i with inhomogeneity parameters Ck+i^ ... ^i, ^ 2 , - - - £fc-i- 


We will also use the following result: 

Lemma 5.2. The trace of the inverse of the quasi-periodic monodromy matrix (15.51) evaluated at some inhomo¬ 
geneity parameter £ n , n € {1,... , N}, is equal to the inverse of the quasi-periodic transfer matrix evaluated at 
fn> i-1- 

trofM^^n)- 1 ] = T {Y \in)~ 1 . (5.7) 
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Proof. For any given numerical matrices Z and Z, we independently compute the two traces tro \Zq -Mo (£«.)] 
and ti'o [Ado(£n) _1 Zq\ as products of ^-matrices or of inverse P-matrices, using the fact that the P-matrix 
(12.11) becomes proportional to the permutation operator P \2 when evaluated at A = 0: Pi2(0|r) = a(0)Pi2- On 
the one hand, we obtain 


tro[Z 0 2Wo(?n)] =tr 0 


k=n +1 


fc -1 


1=1 


Z 0 JJ Rok ( - £fc \T + V X] a l) ‘ «(°) P 0, 

n— 1 / 

xn%(^- €k\T+v X/ a t) j °° 


fc=1 


fc-1 


(=1 


n—1 


k =1 


fc-1 


1=1 


X tl'o 


N , fc-1 

Zo #0fc Kn - £fc \t + V X] a l ) ■ O (°) p 0n 

fc=n+l A 


1=1 

l^n 


= a 


n—l / fc—1 , N / k— 1 

(0) JJPnfcUn - t.k\r + V^af ) T < Z n RnkUn - £k\r + ri^af 

' 7—1 ' T-*— iy~i i 1 ' 7—1 




k=n -\-1 


J=1 

l^n 


where we have used the zero-weight property of the P-matrix, namely [ P n /,.. of, + of\ = 0, as well as the 
commutation relation (12.101) . The second trace can be computed similarly, leading to 


tro[M.«„r I Zo]=^y I \Kl 


k=n -\-1 


fc -1 

tn- £k\r + V^2°l 

1=1 
l^n 


n—l 


Z, T 


ILL 


fc=i 


fc-i 

£n-tk\r + vJ2 a i 

i=i 


hence the result for Z = Z 1 = V. 


□ 


Remark 5.1. It follows from Lemma 15721 and from the inversion relation (12.161) that 

TQefme^) _ e~^ s (?(T) T (K) (d 0) )T (K) (^ 1 1} ) _ e~^ s fl(r) 

A (d 0) )D(d 1} ) ^ T + ? ? S )’ A(ef)D(^ 1} ) <Kt + 1J S)‘ 

Similarly, building on the inversion relation (12.161) and on the proof of Lemma [5721 for different choices of Z , Z, 
one can deduce several other useful identities, such as for instance the following cancellation identities, analog 
to those obtained in If23ll in the non-dynamical case: 

Mj(d 0) ) M ik {& = 0, Vi, j, k = ±. (5.9) 


Remark 5.2. By using similar arguments as in the proof of Lemma 15^21 one can also easily show that 


yr lZM 

M a( &) 


N 

= U{ T - Y f- 


(5.10) 
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Let us now formulate the solution of the quantum inverse problem for the quasi-periodic monodromy matrix 


(1531) . 


Theorem 5.1. The local spin operators E l f (IQ) . understood as operators acting on D^vd).n- can be expressed 
in terms of the entries of the quasi-periodic monodromy matrix (1531) or of its inverse in the following way: 


n —1 


k =n^(6) • • n [T (Y \tk)} ■ v-\ 

k=1 k= 1 

n ^ n—1 

= • [Xi (r) («„)-■]„ ■ n T’d&r" • Tjrt 


k=1 


k =1 


(5.11) 

(5.12) 


Remark 5.3. The reconstructions of Theorem l5 .H are valid on the whole representation space B^vD i.n- on which 
the transfer matrices do not a priori commute. Hence, we have to pay attention to the order in the corresponding 
products. 

Remark 5.4. In 1301 was only formulated the analog of (15.111) in the periodic case. The relation (15.121) is instead 
useful to express local operators in terms of elements of the monodromy matrix with shifted inhomogeneity 
parameters. In fact, using the inversion relation (12.161) for the monodromy matrix, the relation (15.121) can, in the 
periodic or antiperiodic case, be respectively rewritten as 


k = (-i)^ fine*) 


M ^)]- . s 8 (t + V S) 


n— 1 


k= 1 


A«f)D ({«) 

-(-!)¥ nr w (&) 


k =1 


A(^° ) )D(^ ) ) 


9(t) 

) (.Zn ) )]_ i _ j j yriS b(T + pS) 


fc=i 


9(t) 


nin&r -t- 

n- 

■ n ■ T r‘. 


k =1 


(5.13) 

(5.14) 


Proof Let us first show (15.111) for n = 1. The proof is based, as usual (see If26l[34l and the proof of Lemma 1531) . 
on the crucial fact that the /i-matrix (12.11) becomes proportional to the permutation operator P\o when evaluated 
at A = 0. Expressing as in Lcmma [3~2l thc matrix element [A4 fV 1 (^i)]ji as a trace tro [ A/lJ,' (^) E 1 ^ over some 
auxiliary space 0, representing A4 q ; (£i) in terms of a product of ^-matrices and moving in this expression the 
operator E 1 ^ from right to left, using successively that T r°E 1 ^ = E^ TB :i + * 3 , that Pq\E 1 ^ = E\ J Pq\ and that 
(r + pol)E^T l -f J = EfT^P\r + pcrf), we get 

[3f( y ) (6) ] .. = E P ji-i (5.15) 

The general case can be deduced from the case n = 1 by means of Lemma 15.11 Let us express part of the 
product in the left hand side of (15.111) as a multiple trace over auxiliary spaces: 

71—1 

l\T {Y) ^k)-[M^^ n )] ji = tr aia2 ... an 

k= 1 

= t r aia 2 ■■■a r i 


k =1 

n n^ / k—1 

n y «« ii m « u^+i'e< ) ii t b - E i 


_ k=l k =1 


1=1 7 k =1 
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in which we have used (12.101) . The product of monodromy matrices can now be rewritten as in Lemma [5711 and 
reorganized thanks to (12.101) such that 


n— 1 


O r(y) (efc)*[^ (y) (en)], i = t r , 


a\a,2...a n -\ 


k =1 

X 


n— 1 

FL 


a k 


M 


L k =1 

n —1 


(Y) (t \ T + „ v<rf \r Kn E ij 

a n , n n+1... Na\a 2 -.. a„_i l ' ' I / , °l J r “dn 

V (=1 ' 

n—1 

X n M “fc^^+l- Ar “ 1 “2...a fc - 1 (^k + 7?^CTf + T7 <7* ) 
fc=l 


fc-1 


ri—1 \ n—1 


(=1 


(=fc+l ' fc=l 


(5.16) 


The trace over a n can now be explicitly computed similarly as in (15.151) : 

/ n_1 \ 

,{Y) Al r + „^Ajf”El 


tr„ 


M 


a n ,nn+1... Nai < 22 .. .a n —1 


1=1 


— E^ T i_J ' tr 

— rj n 1 T iL an 


M 


00 

a n ,ti n+1... Nai a 2 ■ 


n— 1 


-a*-! Un\T + vJ2 a f) T 


1=1 


(5.17) 


so that the product En T l J can be moved out of the trace from the left in (15.161 ). The remaining multiple trace 
can then be re-expressed as a product of transfer matrices (using again Lemma l5TI) . leading to (15.111) . 

The proof of (15.121) can be performed in a similar way. Considering first the case n = 1 we obtain, as in 
(15.151) . that 


troK-Ad^i)- 1 ] = trofA^&r 1 ] T ^ = T^i)" 1 V T ~ j E?, 


,00/* \-l 


(5.18) 


where we have used Lemma [5721 The general case can be proven by means of Lemma 15. H and Lemma [5721 by 
computing 


n— 1 


[^ (y) (U%rn[ r(y) (^ 1 = t1 ’ 


CLlCL2...CLn 


k =1 




2 X —1 
k 


k =1 


(5.19) 


Reorganizing the product inside the trace and using identity (15.61) for inverse monodromy matrices, we obtain, 
by a similar reasoning as in the previous case, that 


n— 1 


k =1 


k =1 

^ n 

=n 




k =1 


which ends the proof of Theorem l5.ll 


(5.20) 

□ 


As a consequence, one gets 
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Corollary 5.1. The local spin operators Ef + , E n , and the local height operators 5f "'\ 1 < n < N, admit in 
- (q c/n) 

ID(6VD) ^ the following reconstruction in terms of the entries of the n-twisted antiperiodic monodromy matrix: 


n —1 


E n + = II rW fe) • KC ^n) ■ n [^(6)] \ 


( K ), 


k =1 


k =1 


k =1 
n —1 


<«)^ A « # (£n ~ V) 

det q M(£ n ) 


n—1 


n r (&) ■ 


n TA&r 1 , 


- = n ■ n [t i “ j (&)] 


fc=i 
nd. K )ft Nl -1 


fc=l 


fc=l 


n^-OTT-nV’&r 1 ' 


/c=i 

n—1 

2( n ) _ I I t = ( k ) 


det 9 M(£ n ) 


n— 1 




i -i 


fc=i 


fc=i 


(5.21) 

(5.22) 

(5.23) 

(5.24) 

(5.25) 


Proof. The representations (15.21I )- (I5.241) follow directly from (I5.11I) - (I5.12I) taking into account Remark l5T4l and 
the fact that we restrict ourselves to The relation (15.251) follows, as in Theorem 3.2 of | [32l . from a 

trivial recursion on n using the fact that 




(5.26) 


and the solution of the inverse problem for E+f x and E n _ l 


□ 


6 Form factors of local operators 


We are now in position to compute matrix elements of local operators (spin and height) between eigenstates of 
the K-twisted antiperiodic transfer matrix. We obtain the following result: 

Theorem 6.1. The matrix elements of Ef + and Ef~ between generic (T'[| and left and right eigen¬ 
states o/T^ (A) on the following determinant representations: 


(vpW|£++|vpW) = 


wx _ ngTit(^) 


nLit'(6) n+i 


det [5t,t'(?n)] 


flZLi *(&) detjy+i [<St',t(gn ~ rj)] 
nb=it'(6) det 9 M(£ n ) 




u'^rn 

nLi *■'(&) n+i 


det [St',f(U] 


ULi *(&) detjy+i [<%(£n ~ rj)] 
detg M(£ n ) 


( 6 . 1 ) 

(6.2) 


where «Se t'(£n^), e € {0,1}, is an (N + 1) x (N + 1) matrix which corresponds to the matrix Typ (14. 1 6b of the 
scalar product (14.151) with one additional line and one additional column: 


[*M f) )] a , b =[K?] a , b , 

[<% (£i e) )] a , N+1 = Ax,y (£„) qS qp.i, for 

[5t,F(ei e) )] N+1 , 6 = - £o), for 

N+1,N+1 = ^ 


a, b € {1, 

■ ■ • > N}, 

(6.3) 

a € {1,.. 

• j N}, 

(6.4) 

b € {1,... 

j N}, 

(6.5) 



(6.6) 
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The representations (16. lb and (16.21) straightforwardly follow from the solution of the quantum inverse prob¬ 
lem (15.2 1I )- (I5.241) . which enables us to write 


(*?>\E++\*W) = k 


(«)\ _ ^ nt; t (&) , ( K ) 


(^ic(ui 4 k) ) = -« 


rKLi m) 




n6=it'(6) det,M(£ n ) 


- n&=l t(£b) /^(kOi p>(p \ |U/( K )\ — _ K 

' ' 1 ” 1 P> IE=lF&)' ' 1{(n)l r) ~ nCiP(&) 


and from the following lemma: 

Lemma 6.1. The matrix elements of the operators B(£,n' > ) and C(fn'), e € {0,1}, between eigenstates | 
and | Tpd) of the antiperiodic transfer matrix are given by the following determinants: 


<^ K) | B{& |4" } ) = k det fe(£«)], <^1 C(&) |4" } ) = k- 1 dot [*?(&)] , (6.7) 


where <St,t'(£ri^ )> e G {0,1}, is the (N + 1) x (N + 1) matrix with elements (I6.3I) - (I6.6I ). 


Proof From (14.41 ). (13.211) . one can easily compute the action of C(£^), e = 0,1, on a left -eigenstate: 


<*hc& e) )=x; e 


Jy(i 


— £n + £cE) TT ~ & ) 


a=1 he{0,l} N 
h a =0 


n 

b^a 


0(*o,h) 9(^ ha) - 


IM 

x Yl {e iyr)hb K hb qg, b) ) det [0(°’ h )] (0, T+h| . (6.8) 

6=1 N 


Then, using (14.31) and (13.321) . one gets 




r(«) 




,(0) „(1) 


0=1 he{o,i} N 

ha = 0 

N 

b^a - 


Z lyV A x,y (£b) <h b ) (h b ) 

D n b -p) ) 


(-i) N+ “ det le<°;S>(tf>) 


(6.9) 


In this expression, the N x N matrix 0^’'j|j| (fn ^) is obtained from 0( o h) by eliminating the a-th row (containing 


the elements 1 (fa 1 "— £o), 1 < /> < N ) and inserting a new row at position N with elements i iff 1 — £q)> 
1 < 6 < N. Indeed, it follows from (IA.7I) that 


(0 


(-D N+ “ d N et [ e fS] «£>) 




0(*O,t 




l a e(d ha) -al hb) ) N 


( 6 . 10 ) 


It remains to notice that the sum in (16.91) corresponds precisely to the development of the determinant of the 
(N + 1) x (N + 1) matrix Sipiffl) w.r.t. the column N + 1. The proof for the other formula in (16.71 ) is 
similar. □ 
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Theorem 6.2. The matrix elements of local heights operators 5i n fixing the value of the height s £ {/o+r/5; s £ 
{0,1,, N}} at a given site n between generic (\P^| and left and right eigenstates of T^ K \\) on 

- /Q C/1Z) 

®(6VD) N can ^ e wr ^ tten as the following sum of N + 1 determinants: 


("W I*= TT ^ AL_ T det 

' t ' * 1 l ' ' II t'(4) N + 1 ^ n 

6=1 VS0 ' j =0 


( 6 . 11 ) 


where J-~, is the N x N matrix with elements 




n+t ( e in xh 


h =0 


D(£a - V) 


&)• 


( 6 . 12 ) 


Proof Using (15.251) . (14.31 ). (14.41 ) and 13,331) . one obtains 


n_1 m 




Ds Jo,h 


n 




Ax,y (fa) 

D(?a - rj) 


-(NO (ha) 
9t,a 9F,a 


6=1 hG{0,l} N a=l - 

Noticing that to.h £ (h) + tjk', k £ {0,1,..., N}}, one can rewrite 5 Sf t 0 h as 


det r@^ 0,h ^l. (6.13) 

N L J 


J s,t 0h 


1 A 27Ti j(3 r to ’ h) 1 -A - 2 A 2ni]ha 

\ g + N + 1) = - \ g r,(N + l) 1 [ g N+l 

^ N +1 ^ 11 


N + 1 


3=0 


3=0 


a= 1 


which leads to (16.11I )- (I6. 121 ). 


(6.14) 

□ 


To conclude this section, let us briefly comment about these results. Although the matrix elements of local 
spin operators can quite straightforwardly be expressed in terms of a single determinant as in Theorem l6.ll hence 
generalizing the simpler (non-dynamical) six-vertex case | [39l . the situation seems slightly more complicated for 
the matrix elements of local height operators. The latter can nevertheless be expressed as a sum of determinants 
as in Theorem 16.21 The number of terms of this sum being related to the size N of the model, this may (or 
may not, depending on the behavior of the different terms) be a problem for the study of the thermodynamic 
limit. Note however that our study concerns the unrestricted SOS model, for which the heights s arc a priori 
allowed to take an infinite number of valued. In the literature, one usually consider restricted models, such 
as the ABF model f2j] or the CSOS model 1291 1461. for which the crossing parameter rj of the model is rational 
(ij = r/L) and the heights s are only allowed to take a finite number (L) of values. Our present study does 
not directly apply to the ABF case (which corresponds to the case for which there may be poles in (12.21) 1 but 
could easily be adapted to the study of the CSOS case. In the latter case, the matrix elements of local height 
probabilities would be reduced to the sum of only L terms. Although the structure of the determinants at stake 
are a priori quite different, the situation is somehow similar to what happens for the periodic model, which can 
be studied by means of ABA: matrix elements of local spin operators between Bethe eigenstates of the transfer 
matrix can be expressed (at least in the case of the CSOS model) as a single determinant If30l . but it seems that 
matrix elements of local height operators, related to local height probabilities, can only be expressed as sums of 
determinants Il32l . 

12 In fact, due to the antiperiodic boundary conditions that we consider here, the heights of the model are only allowed to take N + l 
values. The consideration of such boundary conditions indeed reduces the actual space of states of the model to the finite-dimensional 
subspace °f the whole representation space N . 
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7 Conclusion 


We have here studied the antiperiodic dynamical 6-vertex model in the SOV framework, for different configura¬ 
tions of the representation space corresponding to the different possible values of the global shift to of the heights 
of the model by half of the periods of the theta function, associated with a couple of parameters (x, y) £ {0, l} 2 . 
We have diagonalized the corresponding antiperiodic transfer matrix, hence obtaining a complete characteriza¬ 
tion of all eigenvalues and eigenstates in terms of a system of discrete equations involving the inhomogeneity 
parameters of the model. We have discussed the rewriting of this characterization in terms of functional equa¬ 
tions of Baxter’s type, and notably in terms of certain classes of solutions of the usual homogeneous function 
T-Q equations. We have also obtained determinant representations for the form factors of the model. 

Several interesting problems remain to be solved. For instance, we have shown the complete equivalence 
between the SOV discrete characterization of the spectrum and eigenstates and the reformulation in terms of 
solutions of the homogeneous T-Q equation (i.e. in terms of Bethe-type equations) in the case of an even 
number of sites only. We plan to consider the case of an odd number of sites in a further study. We also expect 
to be able to use this reformulation so as to consider the homogeneous and thermodynamic limit of the form 
factors formulas that we have obtained here, similarly as what has been recently done in the XXX case If25l . 

Finally, we would like to mention that our results can be used to study, in the SOV framework, the XYZ (or 
eight-vertex) model with various types of quasi-periodic boundary conditions (related to the values of x and y). 
This interesting problem will be the subject of the paper PiTl . 
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A Theta functions, elliptic polynomials and useful identities 

In this paper, 9( A) = 0\ (A|w) denotes the usual theta-function 12011531 with quasi-periods ir and tvlo ( c 3uj > 0), 

OO 

d(z) = -i J (-i)fc e i™(fc+|) a e 2i(fc+|)* j (A.l) 

k=—oo 

OO 

= 2e in % sin Z JJ (1 - e 2 *(«™-*)) _ e 2 i(nnu,+z)j ^ ? ( A . 2 ) 

n =1 

which satisfies 

9(z + tt) = -9(z), 9 (z + ttu) = -e~ ivuj e~ 2iz 9(z). (A.3) 

Throughout the paper, we use the following terminology fi~6l;48l, 44], 

Let r = r^’™) = vrZ + tto-’Z. Let y : F —>• C x be a group homomorphism. We say that a function / of 
z £ C is a theta function of quasi-periods (tv, tvui), of order n and character y if / is a holomorphic function 
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satisfying the quasi-periodicity properties 

/0 + tt) = x(tt) f(z), f(z + tvuj) = x(vrw) e - m ( 2z + 7ra ') /(-). (A.4) 

We say that / is a theta function (or an elliptic polynomial) of quasi-periods ( 7 r, 7ru;), of order n and norm ay 
if / is a theta function of order n and character given by \( 7T ) = (— 1)" and x(ttco') = (—l) n e 2ia f. This is 
equivalent to the fact that there exist constants Ai,..., A n and C with Ai + ... + A n = ay such that 

n 

f(z) = C\[e{z-\ k ). (A.5) 

k =1 


We have the following properties (see for instance fl6ll44ln : 


1. Let 

dim Qn,^ 


be the space of theta functions of quasi-periods (7 r, 7ru ;), of order n E N and of norm a. Then 

= n. 


2. Let /, g e ©iL™'* which coincide at n points x\,..., x n € C: /(xj) = g(xj), 1 < j < n. If x\,... ,x n 
are independent (i.e. if xy — x 3 T and ]E" =1 Xj — a £ T) then f = g. It means that there exists a 
unique theta function (elliptic polynomial) of quasi-periods (tv, tvuj), of order n and of norm a with values 
f(x 1 ),..., f(x n ) at the respective independent points x 1 , x n . It is given by the following interpolation 
formula: 


tm = E 

3 = 1 


9(a ~ EL 1 x fc + xy ~ L TT ~ x k) 


d ( a - EL 1 **) 


n 


^ - *fc) 


f( X j 


(A.6) 


3. Let (L} 1 <j<n be a basis of ©L-Lx Then, for any (x 1 ,..., x n ) G C n , the determinant of the matrix 

(tij(xi))i<i,j,<n is of the form 

n 

t<|e<n ['L( x '*)] = c ■ 9 ( ^X 1 - a] • JJf9(xj - Xj), (A.7) 

J=1 i<3 

where C is some constant. 


We also recall Frobenius determinant formula, for any n-tuples (x] ,..., x n ), (jj\ .... ,y n ) G C n (with 
x* ~~ Vj ^ r, Vi, j) and any t € C (with t £ T): 


det 




0(Xj 

0(Xj 




°{ ELiLi - yj) +1) rii<,<j<n ~ x i) - v*) 

9(t) n::,=i ®(*i - %■) 


(A.8) 


B Inhomogeneous Baxter equation as reformulation of SOV spectrum 

In this appendix we explain how one can show the equivalence of the SOV discrete characterization of the 
spectrum of Theorem 14.11 with the description in terms of elliptic polynomial solutions, with quasi-periods 
(tv, tvcj), of some particular functional T-Q equations with an extra inhomogenous term. 


44 








As explained in Section 14.21 one can modify the functional equation with respect to (14.221) so as to force 
this equation to admit elliptic polynomial solutions with quasi-periods (7r,7ru;). This means introducing some 
gauge transformation as in (14.241 ). so as to adjust the quasi-periodicity properties of the three terms t(A) Q( A), 
a(A) Q( A — ?y) and d(A) Q( A + rj) for t(A) satisfying (14.17l) - (14.18l) . A possible (and somehow minimal) way to 
do it is to choose the function /(A) = fjf\ A) as in (14.26b in terms of two parameters (3 and //. Of course, the 
sum of the aforementioned three terms does not in general cancel, so that one should also add to the equation 
an inhomogeneous term as in (14.251) . As explained in Section l4~2l this still enables ones to recover the condition 
(14.21 ) as long as this inhomogeneous term cancels at all points ^ n \ n £ {1,... , N}, i.e. contains the factor 
A (A) d(A). 

For the class of Q(X) of the form (14.26b with M = N, the function F( A) = F^q{ A) appearing in the 
inhomogeneous equation (14.251 ) is defined in terms of /3, // and Q( A) as 


7 W(X) = P "(-!)" ' * 1 2 3 e 

&(to,o + olq ~ Y,k & + N77) o(p - t 0 , 0) 

| 3e^ x+ ^d(t 0 , 0 ) _ Q(p) fl( A- 

6»(y7ra; - f 0 ,o - + Efc ~ Nr ?) A (F ~ v) 


d + V + yTruj - top - & ~ Nt?) 

< 9 (A — /u + 77) 


(B.l) 


with a.Q = being the norm of the theta function Q( A) of order N. Note that (l4.25b -i lB. 1 1 ) for M = N can 

be seen as an elliptic generalization of the trigonometric inhomogeneous functional equation that was obtained 
in ll42ll . Indeed, under some simple assumptions on the functional dependence of the zeros of the t(A) £ X T , 
when taking the XXZ limit ui —> +ioo (see Remark 12.21 ). one indeed recovers the equation of Theorem 4.1 of 
Il42l . The elliptic analog of Theorem 4.1 of fl42ll can then be formulated as follows: 

Theorem B.l. Let us suppose that the inhomogeneity parameters £ 1 ,..., £n satisfy (13.8b and let us set M = N. 
Then the following two propositions are equivalent: 


1. t(A) is an eigenvalue function of the antiperiodic transfer matrix T(A) (i.e. t(A) £ 

2. t(A) is an entire function of A and, for some f3 £ C \ {0}, there exists a function Q(X) of the form (14.261 ) 
such that (Q(£j), Q(fj — 77 )) 7 ^ ( 0 , 0 ), 1 < j < N, and that t(A) and Q{ A) satisfy the inhomogeneous 
functional equation (14.25I) - (14.26I) - (1B.1I ). 

Ifq £ C \ R, these propositions are also equivalent to: 

3. t(A) is an entire function of A and, for any (3 € C \ {0}, there exists a function Q( A) of the form (14.261) 

such that Q(fj — ??)) 7 ^ ( 0 , 0 ), 1 < j < N, and that t(A) and Q{\) satisfy the inhomogeneous 

functional equation (14.25 1) - (14.26I) - (1B. II ). 


Proof of Theorem 1/1. / I Obviously [7] implies \2\ 

So as to prove that|2] implies[7], let us suppose that, for some 3 £ C\{0}, there exists an elliptic polynomial 
Q( A) of order N such that the function t(A) defined as 


t(A) = 


/?(A) A x , y (A) Q(A - q) + [fjf \A + 77)] 1 d(A) Q(A + 77 ) - a(A) d(A) F™( A) 

Q( A) 


(B.2) 
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is an entire function of A. Then the function e* yA t(A) is a theta function of order N and of norm a.\ = 6 + 

fo,o- Moreover, the particularization of (IB .21) at the 2N points 6 and 6 — //, 1 <j< N, gives 

Q( 6 ) 66 ) - flfHSi) Ax,y ( 6 ) Q (6 - r?) = o, Q (6 - 7,) t (6 - 77 ) - Q(6) = o, (B.3) 

6 ( 6 ) 

for each j G {1,..., N} which, provided that (<2(66 <2(6 —7/)) / (0, 0), means that the matrix (14.191) has 
zero determinant and therefore that t(A) satisfies (14.21) . Hence t(A) G E^=. 

Let us now prove that [7] implies |2] and, in the case 77 ^ M,0 Let d G C \ {0} and let t(A) G For any 
elliptic polynomial Q(A) of degree N and F^q(X) defined in terms of Q(X) by (IB. II ). the function 

(t(A)Q(A) - fjf >(A) Ax,,(A) Q[\ - tj) - D(A) Q(A + v ) + a(A) d(A) (A) 1 (B.4) 

l /A (A+ 7?) J 

is a theta function of order 2N and of norm X^!t=i 6 + hi.o + olq, where oq = A ; is the sum of the roots 

of the elliptic polynomial Q( A). Then the equation (14.251) is satisfied for t and Q if (and only if) it is satisfied in 
N independent points, namely 

• for A = fj, j = 1,... N: t( 6 ) 0(6 -fi + t 0 ,o) <2 ( 6 ) = /3 _ 1 e _iy 6 a x , y ( 6 ) 0(6 - //) <2(6 - 77 ), 

• for A = 6 - 77 , j = 1 ,... N: t (6 - 77 ) 0(6 - p) <2(6 - v) = P ( ' y? 'hf, - p) 0(6 ~ M + 6 , 0 ) <2(6)’ 


provided that 
N 

6 c + 6,0 — Q!Q ^ T. (B.5) 

fc=l 

Since t(A) G £ 7 = satisfies (14.21) . the above system is therefore equivalent to the following system of N equations: 

t(6 ) 0(6 - f + 6,0) <2 (6) = P l( ' yij A x . y (6) <?(6 - m) <2(6 - v), j = 1 ,..., n. (b. 6) 


In general, an elliptic polynomial of order N and of norm olq is completely characterized by its values at N 
independent points. Hence it can be written in the following form: 


o(A )=e ^ ( 66 + f Efi 6 gg) n #64 <?&), 

^ 0(26 6-« q ) 


k =1 


e=i ~ 6 ) 

i+k 


(B.7) 


provided Yht 6 — °Q £ r. Hence, the system (IB.61 ) is in fact a system of N homogeneous linear equations in 
the N unknowns Q( 6 ), n G {1,..., N}, which can be written as: 


N 

^[C- t (!3,a Q )\ jk Q(t k ) = 0, 


k =1 


j = 1,... , N, 


(B. 8 ) 


where Ci(j3, olq ) is the N x N matrix of elements 


\Ct(P,a Q )\ ab = 5 a b 


P e^ a t (6) 0(6 ~ M + t 0 ,o) 0(6 - & - »7 + £/ & “ «o) jj 0(6 ~ 6 - v) (B 9) 


^x.y (6 


0(6 - aO 


0(XA6 - «q) 


1 =1 
^6 


0(6 - 6 ) 
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Note that, using Frobenius determinant formula (1A.8I) and the formula for the determinant of the sum of two 
matrices, one can express the determinant of the matrix ( IB.91 ) in the form 


N 

det [CtOMo)] =£(-ir/? N - 

n =0 


~ a Q~ n v) 

0(Ye& - «<?) 


x 


e n{ 

PC{1.N}a£P L 

#P=n 


e *y$a F + fp.o) 

Ax,y(£a) ~ p) 


n 

b£P 


0(ga ~ & + V) \ 

~ &>) j 


(B.10) 


where the second summation in (IB. 101) runs over all subsets P of the set {1,..., N} with cardinality n. 

The system (IB. 81) admits some non-zero solution (<2(£i),..., Q(£n)) if and only if the determinant of the 
matrix ( IB. 91 ) is zero. At P = 0, this determinant simplifies into 


det [Ct(0, &q)] 
INI L J 


/ pN ®( a Q ~ Ye& + N, l) 
e ( a Q~Ye&) 


so that 


(B.l 1) 


det [Ct(0 , Q!q)] = 

ddet N [Cj(P,a Q )\ 
da Q 


0 3(fci, /C 2 ) G Z 2 , a.Q = ^ 

£ 


QQ=J2i £e—Nr)+Trki+Truk2 
f3=0 


7^0. 


& — Nr) + irki + nujk2, 


(B.12) 

(B.13) 


Flence we can apply the implicit function theorem for holomorphic functions: V(A: 1 , ko ) € Z 2 , there exist some 
open vicinities and V(k 1} k 2 ) °f 0 and of a|°| k ^ = Ye& ~ Np + irki + Truk 2 respectively, and there 

exists a unique holomorphic function a(ki,k 2 ) '■ ^(ki,k 2 ) Y k ito) a (fci,fc 2 )( 0 ) = k 2 ) suc ^ that 


{(/3,«q) G U(k\,k 2 ) ^ V(k\,k 2 ) I det [Ct(/3, Q:q)J Oj’ ^ {Pi ^(k\ ,k 2 ) (/^)) I P G i3(fci,fc2) j" ■ 


(B.14) 


Hence the system 


N 

^[C- t {l3,a {klM) (P))] jk q k = 0, j = 1,..., N, (B.15) 

k =1 

admits, for all f3 G U( k] k2 y a non-zero solution (<[\..... q^). Due to the form of (IB.151 ). one can choose this 
solution such that all qj = q-j(8) are continuous function of (3 in Ur kl k2 y 

Note that, at (3 = 0, we have qj( 0) qP 0 for all j G {1,..., N} (it is clear from the system (IB.61 ) that 
all the roots of the solution Q{ A) are, up to T-periodicity, at the points — r], 1 < j < N, when 8 = 
0). Hence it is always possible to choose U( klk2 \ such that qj(/3) / 0 for all j G { 1,... , N} and for all 
P in U( klM ). Moreover, since a^ klM) {0) - Ye& t T and a^ klM) {0 ) ~Ye&~ h,o i T, and since the 
function a( fcl fe2 )(/d) is holomorphic, one can also choose U( k] k2) such that o-( klt k 2 ){P) ~ Ye & ^ I an d that 
a(fci,fc 2 )(/ ? ) -Ytil- *o,o i r for any P G Uy. lM) . Hence we have showng] 

Let us now notice that, up to a similarity transformation, the matrix Ci(P, qq + mu) is proportional to the 
matrix Ci(pe~"8 olq). Hence the system (IB. 81) admits a non-zero solution for (p. olq) if and only if it is the case 
for the system corresponding to (/3e 2m7? , o,q + niruj) (Vn G Z). 
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Let P £ C \ {0}, and let us suppose moreover that rj ^ R. Hence, there exists n £ Z such that (3 e 2lkri £ 
U m , so that the system associated with (/3e 2m?7 , cq 0 , 0 )(/3e 2mr? )) admits a non-zero solution such that each 
unknown (i.e. each Q(fj) solution to this system) is itself non-zero. It follows from the previous remark that 
the system for (/3, a( 0 ,o) (/3e 2mr? ) — riTTUj'j admits also a non-zero solution which is such that Qifj) f 1 0 for each 
j £ {1,..., N}. Hence we have shown|3 □ 

The inhomogeneous functional equation (14.251) - (14.26I) -( 1B . 1 1) for M = N is not the only functional equation 
that can be considered in this framework. For instance, it is also possible to completely characterize the SOV 
spectrum in terms of a functional equation of the same type (still using (14.261 ) to define the gauge function 
/(A) in terms of some arbitrary parameter p) but for an elliptic polynomial Q( A) of degree M = N + 1 and 
arbitrary norm gq . Of course the inhomogeneous term has to be adapted accordingly and it appears slightly 
more complicated in this case. In that way, increasing the degree of the elliptic polynomial Q(\) corresponds 
to increasing the number of free parameters in the inhomogeneous equation: apart from the gauge parameter 
P, we have one free parameter (the parameter p) for the degree N, and two free parameters (p and gq) for the 
degree N + 1. Hence, the choice of the degree N — 1 for Q{ A) seems to be the minimal possible if one considers 
equations of the form (I4.25I )- (I4.26I ): in that case we have still N unknown parameters which are the parameter 
p appearing in the definition (14.261) of /(A), as well as the N — 1 roots of Q{ A), and there does not remain any 
free parameter (except /3), cf. footnote [ 6 ] 

It is interesting to remark that, if we have a complete description of the transfer matrix spectrum in terms 
of the elliptic polynomial solutions Q(X) of degree M of some inhomogeneous functional equation of the form 
(14.251 ) with associated function /(A), it is possible to rewrite the transfer matrix eigenvectors in a generalized 
Bethe form, in terms of the roots of the corresponding elliptic polynomial Q(X). More precisely, defining the 
states, 


i</> 

(^Mjl 


N 

e n 

he{o,i} N a=l 


e iyr? A x , y (e a ) 

KD(£ a - rf) 


/(&)} dct[0(°- h )] |h,-M), 


N 

e n /&*)) la d N et [ 0(o,h) ] < m > h 

he{o,i} N a=l 


(B.16) 

(B.17) 


we have the following result: 


Corollary B.l. Let t(A) be an eigenvalue of the K-twisted antiperiodic transfer matrix (i.e. t(A) F E.p), and let 
M 

Q(A) = n 0 (A“Aj) (B.18) 

i=i 

be such that t(A) and Q( A) satisfy the functional equation (14.251) for some function /(A). Then the T( A)- left 
and right eigenstates with eigenvalue t(A) can be represented as 

M M 

l^ K) > = n [ e * yT ^( r )- 1 ^ , (Aa)] <^1 = (</l n [^(A(B.19) 

a= 1 ci— 1 

where the order of the operators in each bracket [...] has to be kept as it appears. 


In other words, if /(A) is fixed, which is the case for the complete characterization of the transfer spectrum 
that we have obtained in Theorem lB.ll in terms of elliptic polynomials solutions of (14.25l) - (14.26l) for M = N, then 
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(IB . 1 61) and (IB. 171) have to be understood as some fixed pseudo-vacuum states. The corresponding eigenstates 
are then obtained by multiple action, on these pseudo-vacuum states, of the (slightly dressed) operator T>( A) 
evaluated at the roots of the Bethe equations, i.e. in a form very similar to what happens in the context of ABA. 

Let us finally insist on the fact that the possibility to write such an ABA-type representation for the SOV 
transfer matrix eigenstates is very general for models solved by SOV, as soon as we have some characterization 
of the transfer matrix spectrum in terms of (generalized) polynomial solutions to some (homogeneous or inho¬ 
mogeneous) T-Q functional equation. It is not restricted to reformulations of the SOV spectrum by inhomoge¬ 
neous Baxter equations; indeed, we can construct similar reformulations, using multiple action of some slightly 
different operator, in the framework of the characterization in terms of the solutions of the homogeneous T-Q 
equation as obtained in Section [4] It is not restricted to the model under consideration; indeed, a similar rewrit¬ 
ing is possible for all the integrable models so far solved in the SOV framework lf43l[37l 12T. L4 liT3lf39ll40l!42ll 
as soon as we have a characterization of the SOV spectrum by a finite system of (generalized) Bethe equations. 
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